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ABSTRACT 



Quantum algebras are a mathematical tool which provides us with a class of sym- 
metries wider than that of Lie algebras, which are contained in the former as a 
^ ' special case. After a self-contained introduction to the necessary mathematical 

. tools (g- numbers, g-analysis, q'-oscillators, q'-algebras), the sug(2) rotator model 

and its extensions, the construction of deformed exactly soluble models (Interact- 
ing Boson Model, Moszkowski model), the use of deformed bosons in the description 
of pairing correlations, and the symmetries of the anisotropic quantum harmonic 
1/^ , oscillator with rational ratios of frequencies, which underly the structure of superde- 

0^ ' formed and hyperdeformed nuclei, are discussed in some detail. A brief description 

of similar applications to molecular structure and an outlook are also given. 

^ ' 1. Introduction 

Quantum algebras Sffll'l (also called quantum groups) are deformed versions of 
, . ■ the usual Lie algebras, to which they reduce when the deformation parameter g is set 

! equal to unity. From the mathematical point of view they are Hopf algebras H. Their 

\ use in physics became popular with the introduction of the g-deformed harmonic 

oscillator (sec. 10) as a tool for providing a boson realization of the quantum algebra 
sUg(2) (sec. 14), although similar mathematical structures had already been known 
(sec. 11). Initially used for solving the quantum Yang-Baxter equation 0, quantum 
algebras have subsequently found applications in several branches of physics, as, for 
example, in the description of spin chains, squeezed states, rotational and vibrational 
nuclear and molecular spectra, and in conformal field theories. By now several kinds 
of generalized deformed oscillators (sec. 12) and generalized deformed su(2) algebras 
(sec. 17) have been introduced. 

It is clear that quantum algebras provide us with a class of symmetries which 
is richer than the class of Lie symmetries, which are contained in the former as a 
special case. It is therefore conceivable that quantum algebras can turn out to be 
appropriate for describing symmetries of physical systems which are outside the realm 
of Lie algebras. 



Here we shall confine ourselves to applications of quantum algebras in nuclear 

structure physics. The structure of this review is as follows: In order to make this 
review self-contained, we are going first to give a brief account of the necessary tools: 
g-numbers and g-analysis (sees 2-9), g-deformed oscillators (sees 10-13), g-deformed 
su(2) algebras (sees 14-18). The remainder will be devoted to specific applications 
in nuclear structure problems, starting with phenomenology and advancing towards 
more microscopic subjects. The sUg(2) rotator model (sees 19-23) and its extensions 
(sees 24-26), as well as the formulation of deformed exactly soluble models (Interact- 
ing Boson Model (sees 27-29), Moszkowski model (sec. 30)) will be covered in some 
detail. Subsequently, the use of quantum algebraic techniques for the description of 
pairing correlations in nuclei (sees 31-33), as well as the symmetries of the anisotropic 
quantum harmonic oscillators with rational ratios of frequencies (sec. 34) will also be 
considered in some detail. The latter are of current interest in connection with the 
symmetries underlying superdeformed and hyperdeformed nuclear bands (sec. 34). 
Finally, a brief account of applications of the same techniques to molecular structure 
(sec. 35) and an outlook (sec. 36) will be given. 



2. g- numbers 

The q-number corresponding to the ordinary number x is defined as 

N = ^ff, (2.1) 

where g is a parameter. The same definition holds if x is an operator. We remark 
that g- numbers remain invariant under the substitution q — > q~^. 
If q is real, g-numbers can easily be put in the form 

[x] = (2.2) 
smh(r) 

where q — e^ and r is real. 

If g is a phase factor, g-numbers can be written as 

, , sinfra;) 

N = ^Vt' 2.3 
sm(r) 

where g = e*"^ and r is real. 

In both cases it is clear that in the limit g — > 1 (or, equivalently, r — > 0) g-numbers 
(or operators) tend to the ordinary numbers (or operators): 

lini[a;] = x. (2.4) 
A few examples of g-numbers are given here: 



[0] = 0, [1] = 1, [2] = g + g-i, [3] = g^ + 1 + g-^. (2.5) 



Identities between g'-numbers exist. They are, however, different from the famihar 
identities between usual numbers. As an exercise one can show (using the definition 
of g-numbers) that 

[a][b + 1] - [b][a + 1] = [a - b]. 
The q-factorial of an integer n is defined as 



N! = [n][n-l]...[2][l]. 
The q-binomial coefRcients are defined as 



[m]\ 



[m — n\\ [n] ! ' 
while the q-binomial expansion is given by 
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In the limit g — > 1 we obviously have 

[n]\ —>■ n\ and 
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(2.6) 
(2.7) 



(2.8) 



(2.9) 



(2.10) 



where n\ and are the standard factorial and binomial coefficients respectively. 

It should be noticed that two-paremeter deformed numbers have also been intro- 
duced 

[xU - ^r^- (2.11) 
In the special case p = q they reduce to the usual g-numbers. 



3. gr-deformed elementeiry functions 

In addition to g-deformed numbers and operators, g-deformed elementary func- 
tions can be introduced. The q-exponential function is defined as 

oo n 

= E OT^"' (3-1) 

n=0 l-'^J- 

while the q-trigonometric functions are defined as 

oo ™2n+l 
oo 2n 

cos,{x) = fyrj^- (3-3) 



It should als 

polynomials 

34.3). 



; noticed that g-deformed polynomials, Buch as g-deformed Hermite 
and g-deformed Laguerre polynomials also exist (see also subsec. 



4. g-derivatives 

Proceeding along this path one can build a new differential calculus, based on 
g-deformed quantities (see 00 for concise expositions). For this purpose the g- 
dehvative is defined as 

The similarity between the present definition and the one of g-numbers (eq.(2.1)) is 
clear. 

Using the definition of the g-derivative one can easily see that 

D^(ax") = a[n]x"-\ (4.2) 

Dleq{ax) = aeq{ax). (4.3) 
From the definition of the g-derivative one can further derive the sum rule 

Dlifix) + g{x)) = Dlfix) + Dlg{x), (4.4) 

as well as the rule 

D^axi ± 6x2]™ = ±[m]b[axi ± 6x2]™"\ (4.5) 

where a and b are costants and [axi ±6x2]"^ is given by the g-binomial expansion (eq. 
(2.9)). One can also prove the q-integration by parts formula 

DUfi.Hx)) = (4.6) 

[q — q )x 

From this, the following two forms of the Leibnitz rule can be derived 

Dl{f{x)g{x)) = {Dlf{x))g{q-'x) + f{qx){Dlg{x)), (4.7) 

Dl{f{x)g{x)) = {Dlg{x))f{q-'x) + g{qx){Dlf{x)). (4.8) 
In addition one can show the property 

Dlf{qx)=qDlf{x)U=q,, (4.9) 

and the chain rules 

DlJix) = Inifix) , (4.10) 

Dlfix-) = [n]x-'Di:f{x-), (4.11) 



where a is a constant. Another useful resuh is 



1 n— 1 

Dffix) = - E Dlfig'^^-^^-'^x). (4.12) 



5. g-integration 

The g-integration (see 0'0 for concise expositions) in the interval [0, a] is defined 

by 

/ /(x)rf,a: = a(g-i-g)5:g2"+V(g'"+M, (5-1) 
while for the interval [0, oo) one has 

/ f{x)d,x = {q-'-q) Y: (5.2) 

•^0 n=-oo 

The indeGnite q-integral is defined as 

/oo 

/(a;)rfga; = (g"^ - g) E g^"+^a;/(g2"+^s) + constant, (5.3) 

n=0 

where < g < 1. For entire functions f{x) one can easily see that this g-integral 
approaches the Riemann integral as g — > 1, and also that the operators of g-different- 
iation and g-integration are inverse to each other 

Dl J f{x)d,x = fix) = J Dlf{x)d,x. (5.4) 

One can also easily see that 

ax^~^dqX = -r^ax^ + constant, (5.5) 

[n] 

eq{ax)dqX = -eq{ax) + constant. (5.6) 
From (4.6) one can also prove the following formulae of integration by parts 

^ fiqx)iDlgix))dqX = fix)gix)\:=^, - £ {Dlf{x))g{q-'x)dqX, (5.7) 

f{q-'x){Dlg{x))dqX = f{x)g{x)\:z^o ' P {Dlf{x))g{qx)dqX. (5.8) 
Jo 

The following formulae can also be proved 

f{x)daqX = a f{x)dqX, (5.9) 



/(x")c/qnx" = [n] J x''-^f{x'')dqx, (5.10) 

/ fix)d,x = — 5: g^'^-^"-!) / fiq"'~^--'^x)d,nx. (5.11) 

The q-analogue for Euler's formula for the function r(x) is 
rC 



e,(-x)x"c;gX = M[n - l][n - 2] . . . [1] = [n]\. (5.12) 
A proof of this formula can be found in 0. 
6. Q-numbers 

The definitioiLiif g-numbers given in sec. 2 is not the only possible one. Much 
literature exists llBOJ using the definition of Q-numbers 

H« = (6.1) 

where x can be a number or an operator and Q is a deformation parameter. Q is a 
real number {Q ^ 0, 1). The notation Q = e^, where T a real number, will be often 
used. The subscript Q will be used in this review in order to distinguish deformed 
numbers defined as in eq. (6.1) from these defined by eq. (2.1). It is clear that in 
the limit Q ^ 1 (or, equivalently, T ^ 0) Q-numbers become ordinary numbers, i.e. 

[x]q X. 

A few examples of Q-numbers are given here: 

[0]q = 0, [1]o = 1, [2]Q = g + l, [3]q = Q^ + Q + 1. (6.2) 

Q-numbers clearly do not remain invariant under the substitution Q Q^^. One 
can easily prove that 

Wq = Q^-'Wi/q- (6.3) 
Q-numbers are connected to g-numbers through the relation 

M = ?''"Nq, with Q = g2. (6.4) 

The definitions of Q-factorials and Q-binomial coefficients still look like the ones 
given in eqs. (2.7)-(2.8): 



W- = [n]Q[n-l]Q...[l]Q, 



(6.5) 



n 



As it can be easily seen from eq. (6.1) under the substitution Q 
and 

Qk{n—k) 



1/Q 
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Q-factorials are connected to g-factorials by 



\n 



]! = g^"("-')/'[n]Q!, with Q = q' 



Q ^ one obtains 
(6.7) 

(6.8) 
(6.9) 



7. Q-deformed elementary functions 

The definitions of Q-deformed elementary functions look similar to these given 
in sec. 3. The Q-deformed exponential function is defined as 



ax) 



E 

n=0 



tX 



n 



and satisfies the property 



x)ei/Q{-x) = 1. 



(Notice that eQ{x)eQ{—x) 7^ 1.) 

The Q-deformed trigonometric functions are defined as 



X 



2n+l 



smQ[x) 



[2n+l]QV 



X 



2n 



COSq(x) = ^(-1) 

n=0 ^^"'Jt 



One can easily show that 



sinQ(x) = — {egiix) - eqi-ix)) , 

cosq{x) = ^ {eqiix) + egi-ix)) . 
Instead of the familiar identity sin^(a;) + cos^(a;) = 1 one has 

sing (a;) sini/Q(x) + cosQ(a;) cosi/Q(a;) = 1. 



(7.1) 
(7.2) 



(7.3) 
(7.4) 

(7.5) 
(7.6) 

(7.7) 



Q-deformed polynomials, such as Q-deformed Hermite polynomials and Q-defo- 
rmed Laguerre polynomials also exist Ej. 



8. Q-derivative 

Given the function f{x) one defines its Q-derivative 111 by the expression 

The similarity between this definition and the one of Q- numbers (eq. (6.1)) is clear. 
One can easily prove that 



which looks exactly like eq. (4.2). In addition one has 

D'^eqlax) = aeqlax), (8.3) 

D^(^i/Q{ax) = aei/Q{aQx), (8.4) 

Dosing (ax) = acosQ(ax), (8.5) 

Z)^ cosQ(ax) = — a sinQ(ax). (8.6) 

One can also easily see that sing (ax) and cosQ(aa;) are the linearly independent 
solutions of the Q-differential equation 

{D^Yu{x) + a\{x) = 0, (8.7) 

while the functions sini/Q(ax) and cosi/Q(ax) satisfy the equation 

{D^)\{x) + a\{Q^x) = 0. (8.8) 

The following Leibnitz rules can also be shown: 

{fi{x)f2{x)) = {d^Mx)) UQx) + h{x) {D^f2{x)) , (8.9) 

{h{x)Mx)) = (d^Mx)) Ux) + MQx) {D^f2{x)) . (8.10) 
One can further obtain 



r^oMx) _ jP^M-)) Mx) - h{x) (d^Mx)^ 

^ f2{x) f2{Qx)f2{x) ■ ^ - ' 

For the second derivative of f{x) one has 

(D^Yfix) = {Q- l)-'Q-'x-' [f{Q'x) - (Q + l)f{Qx) + Qf{x)} , (8.12) 



and by mathematical induction we obtain the general formula 



{D^Tfix) = {Q- l)-"Q-"("-i)/22;-« J2 



k=0 



n 
k 



(_l)fcgfc(fc-i)/2j(gn-fca;). (8.13) 



9. Q-integration 

In a way analogous to thatof sec. 5 the definite Q-integral of the function f{x) 
in the interval [0, 1] is defined liHl as follows 

/ f{x)dQX={l-Q)Y.fmQ% (9.1) 



assuming that Q is real and \Q\ < 1, while for the definite integral of f{x) in the 
interval [0, 00], we have 

/•OO ^ 

/ f{x)dQX = {i-Q) Q'lm- (9.2) 

For the indefinite Q- integral of f{x) one has 

/OO 

f{x)dQX = (1 - Q)x Q'fiQ'x) + const. (9.3) 
s=-oo 

One can easily check that Q-differentiation and Q-integration are operations inverse 
to each other 

D« / f{x)dQX = fix). (9.4) 



The formula for Q-integration by parts reads 

;D«/i(x)) f2{x)dQX = h{x)h{x) - I h{Qx) (/^«/2(x)) dQX. (9.5) 



10. The g-deformed harmonic oscillator 

The interest for possible applications of quantum algebras in physics Jiaa^iieen 
triggered in 1989 by the introduction of theo-deformed harmonic oscillator EJ'EJo, of 



which earlier equivalent versions existec 

The g-deformed harmonic oscillator llSll3ll§Ej'El'Bil is defined in terms of the creation 
and annihilation operators and a and the number operator A^, which satisfy the 
commutation relations 

[N,a)] = a\ [N,a] = -a, (10.1) 
aa^ - q^^a^a = q^^ . (10.2) 



In addition the following conditions of hermitian conjugation hold 

(a^)^ = a, = N. (10.3) 

Eq. (10.1) is the same as in ordinary quantum mechanics, while eq. (10.2) is modified 
by the presence of the deformation parameter q. For g — > 1 it is clear that eq. (10.2) 
goes to the usual boson commutation relation [a, a^] = 1. An immediate consequence 
of (10.2) is that 

a^a^[N], aa^^[N+l]. (10.4) 

Thus the number operator is not equal to a^a, as in the ordinary case. The 
operators and a are referred to as q-deformed boson creation and annihilation 
operators respectively. 

The basis of the Fock space is defined by repeated action of the creation operator 
on the vacuum state, which is annihilated by a: 

a|0) = 0, |n>=^^|0>. (10.5) 



The action of the operators on the basis is given by 

N\n>^n\n>, (10.6) 



aV \/[n+l]\n >, (10.7) 



a\n>^ ^[n]\n-l> . (10.8) 

We remark that these equations look very similar to the ones of the ordinary case, 
the only difference being that q'-numbers appear under the square roots instead of 
usual numbers. 

The Hamiltonian of the g-deformed harmonic oscillator is 

H ^—{aa^ + a^a), (10.9) 

and its eigenvalues in the basis given above are 

E(n)^^{[n] + [n + l]). (10.10) 

One can easily see that for q real the energy eigenvalues increase more rapidly than 
the ordinary case, in which the spectrum is equidistant, i.e. the spectrum gets "ex- 
panded". In contrast, for q being a phase factor (g = e'"^ with r real) the eigenvalues 
of the energy increase less rapidly than the ordinary (equidistant) case, i.e. the spec- 
trum is "compressed" . In particular, for q real {q — e'^) the eigenvalues can be written 
as 

fi^j sinh (r (n + ^)) 



while for q being a phase factor (g = e*"^) one has 



Tn^ sin ( r ( n + i , 
E{n) = — ^ ^ (10.12) 



2 sin 2 



In both cases in the limit q ^ 1 (r ^ 0) the ordinary expression 



E{n) = huo[n + ]-] (10.13) 



is recovered. 

In addition, the following commutation relation holds 

[a,a)] = [N +l]-[N]. (10.14) 
For q being a phase factor, this commutation relation takes the form 

(27V+l)r 

[a, at] = I . (10.15) 

cos I 

It is useful to notice that the g-deformed boson operators and a can be expressed 
in terms of usual boson operators and a (satisfying [a, a^] = 1 and N = a^a) 
through the relations ElU 



„ = . E±il„ = JM. „t.„, E±il = .M„t. (10.16) 

VA^ + l VA^ ^ N + 1 ^ N ^ ' 

The square root factors in the last equation have been called q-deforming functionals. 

For q being a primitive root of unity, i.e. q = e^'^'^l^ (fc = 2, 3, . . . ), it is clear the 
the representation of eqs (10.5)-(10.8) becomes finite-dimensional and has dimension 
/c, since only the vectors |0 >, |1 >, , |/c — 1 > can be present. This case has been 
related to the system of two anyons In what follows we are going to assume that 
q is not a primitive root of unity. 



11. The Q-deformed harmonic oscillator 

A different version of the deformed harmonic oscillator can be obtained by defining 
the operators 6, through the equations 

a = q^l%q-^l\ o) = q^'\-^'%l (11.1) 

Eqs. (10.1) and (10.2) then give 

[N,h^]=h\ [AT, 6] = -6, (11.2) 



This oscillator liasJDeen first introduced by Arik and Coon and later considered 
also by Kuryshkin EJ. One then easily finds that 

b^b=[N]Q, bb^ = [N + l]Q, (11.4) 

where Q = and Q-numbers are defined in (6.1). The basis is defined by 

b\0) = 0, \n>=^^=\0>, (11.5) 

while the action of the operators on the basis is given by 

N\n>=n\n>, (11.6) 

b^\n>= ^J[n + l]Q\n + l>, (11.7) 



b\n >= .J[n]Q\n - 1 > . (11.8) 
The Hamiltonian of the corresponding deformed harmonic oscillator has the form 



H = ^{bb^ + b^b), (11.9) 



the eigenvalues of which are 



E(n) = ^i[n]Q + [n + l]Q). (11.10) 

One can easily see that for Q = , where T > and real, the spectrum increases 
more rapidly than the ordinary (equidistant) spectrum, while for Q = e^, with T < 
and real, the spectrum is increasing less rapidly than the ordinary (equidistant) case. 
From the above relations, it is clear that the following commutation relation holds 

[&,6l = g^. (11.11) 



12. The generalized deformed oscillator 

In addition to the oscillators described in the last two sections, many kinds of 
deformed oscillators have been introduced in the literature (see il for a list). All 
of them can be accommodated within the common mathematical framework of the 
generalized deformed oscillator which is defined as the algebra generated by the 
operators {l,a,a^, A^} and the structure function $(x), satisfying the relations 

[a,N] = a, [a\N] = -a\ (12.1) 



a^a = <^{N) = [N], aa'^ = <I>(A^ + 1) = [A^ + 1], 



(12.2) 



where is a positive analytic function with $(0) = and N is the number 

operator. Prom eq. (12.2) we conclude that 

Ar = $-i(ata), (12.3) 

and that the following commutation and anticommutation relations are obviously 
satisfied: 

[a,a^] = [N + 1]-[N], {a,a^} = [N + 1] + [N]. (12.4) 

The structure function $(a;) is characteristic to the deformation scheme. In Table 
1 the structure functions corresponding to different deformed oscillators are given. 
They will be further discussed at the end of this section. 

It can be proved that the generalized deformed algebras possess a Fock space of 
eigenvectors |0 >, |1 >,..., |n >, ... of the number operator A^" 

N\n >— n\n >, < n\m >— 5nm, (12.5) 

if the vacuum state |0 > satisfies the following relation: 

a|0>=0. (12.6) 

These eigenvectors are generated by the formula: 

where 

n n 

fe=i fe=i 

The generators and a are the creation and annihilation operators of this deformed 
oscillator algebra: 



a|n >= y [n]a|n - 1 >, a}\n >= ^[n + l\a\n + I > . (12-9) 
These eigenvectors are also eigenvectors of the energy operator 

H ^ —{aa^ + a^a), (12.10) 

corresponding to the eigenvalues 

E{n) = ^($(n) + $(n + 1)) = ^([n] + [n + 1]). (12.11) 

For 

$(n) = n (12.12) 



Table 1: Structure functions of special deformation schemes 



$(a;) Reference 

i X harmonic oscillator, bosonic 

algebra 

ii q'^-q^^ g-deformed harmonic oscil- 
' ' latorffl 

iii ^fj" Arik-Coon, Kuryshkin^___or 



V 



xn 
xiii 



Q-deformed oscillator 



iv ^ 2-parameter deformed oscil- 

' " lator §§§ 



x{p + 1 — x) parafermionic oscillator S 



vi sinh(rx)jnhHp+i-x)) g-deformcd _ parafermionic 



oscillator 

vii a;cos^(7rx/2) + {x + p — 1) sin^(7ra;/2) parabosonic oscillator 

sinh(Ta;) cosh(T(a;+2Aro-l)) qq^.'^ (^-j^ x / 2) 

+ sinh(r) cosh(r) l^TX/ ZJ dllator H 

ix sin^ 7rx/2 fermionic algebra 

X sin^ 7rx/2 g-deformed fermionic alge- 



viii '''fll}r(x+2NT-K-^cosh(rx) .'2/' g-dcforined parabosoulc OS 



xi ^ generalized g^eformed 



bra iSiiiS'Eil 

generalized 
fermionic algebra 



sn{T) 



one obtains the results for the ordinary harmonic oscillator. For 

^^^^ ^ qn_q n ^ 

q-q 

one has the results for the g-deformed harmonic oscillator, while the choice 

$H = = Mq (12.14) 

leads to the results of the Q-deformed harmonic oscillator. Many more cases are 
shown in Table 1, on which the following comments apply: 

i) Two-parameter deformed oscillators have been introduced HSill, in analogy to 
the one-parameter deformed oscillators. 



ii) Parafermionic oscillators eHI of order p represent particles of which 
number which can occupy the same state is p. Parabosonic oscillators 
introduced. 

iii) g-deformecLyersions of the parafermionic and parabosonic oscillators have also 
been introduced Eaa. 

formed versions of the fermionic algebra have also been introduced 
as well as g-deformed versions of generalized g-deformed fermionic al- 
Lhowever, that g-deformed fermions are fully equivalent 



gebras 0. It has been prcw 
to the ordinary fermions 



e maximum 
can also be 



13. The physical content of deformed harmonic oscillators 

In order to get a feeling about the physical content of the various deformed har- 
monic oscillators it is instructive to construct potentials giving spectra similar to 
these of the oscillators. 



13.1. Classical potentials equivalent to the q-oscillator 

Let us consider the g-deformed harmonic oscillator first. For small values of r 
one can take Taylor expansions of the functions appearing there and thus find an 
expansion of the g-number [n] of eq. (2.1) in terms of powers of r^. The final result 
is 

^2 4 ^6 

\n\ =n±—(n- n^) + (7n~ lOn^ + 3n^) ± (31n -49n^ + 21n^ - 3n^) + . . . , 

^ ^ 6 ^ ^ 360^ ^ 15120^ ^ 

(13.1) 

where the upper (lower) sign corresponds to g being a phase (real) . Using this expan- 
sion the energy of the g-deformed harmonic oscillator of eq. (10.10) can be rewritten 

as 2 2 

E(n)/(M = {n + ^)il±^)T jin + ^)' + • • • (13.2) 
On the other hand, one can consider the potential 

V{x) = Vo + kx^ + Xx^ + fxx^ + ^x^ + .... (13.3) 

If A, /i, ^ are much smaller than k, one can consider this potential as a harmonic 
oscillator potential plus some perturbations and calculate the corresponding spectrum 
through the use of perturbation theory e3 (see also subsec. 32.2). In order to keep the 
subsequent formulae simple, we measure x in units of (^/(2mc(j))^/^. Using standard 
first order perturbation theory one finds that the corresponding spectrum up to the 
order considered is 

Ein) = Eo+(2K+25/i)(n+^) + (6A+2450(n+^)2+20/i(n+^)=^ + 70e(n+^)^ (13.4) 



By equating coefficients of the same powers of {n + |) in eqs. (13.2) and (13.4), we 
can determine the coefficients appearing in the expansion of the potential given in eq. 
(13.3). The final result for the potential, up to the order considered, is 

V{x) = (i ± j)x' T {^x^ (13.5) 

We see therefore that to lowest order one can think of the g-oscillator with a small 
value of the parameter r as a harmonic oscillator perturbed by a term. It is clear 
that if we go to higher order, the next term to appear will be proportional to r^x^°. 

The results of this subsection are corroborated by an independent study of the 
relation between the g-deformed harmonic oscillator and the ordinary anharmonic 
oscillator with x^ anharmonicities eZI. 

13.2. Classical potentials equivalent to the Q-deformed oscillator 

Similar considerations can be made for the Q-oscillator. Defining Q = it is 
instructive to construct the expansion of the Q-number of eq. (6.1) in powers of T. 
Assuming that T is small and taking Taylor expansions in eq. (6.1) one finally has 

= n+'^{n^ -n) + ^C^n^ - + 1) + ^{n^ - 2n^ + n^) + . . . (13.6) 

Then the corresponding expansion of the energy levels of the oscillator of eq. (11.10) 
is 

Ji rp2 rri3 -, rri rp2 f-TiS 1 

Ein)/{huj) =E' + {!-- + + . . + _) + ( + 1 . . .)(^ + _)2 

V ;/v ; ^ 2 8 16 2' ^ 2 A 48 2' 

rp2 nnS I 1 

Comparing this expansion to eq. (13.4) and equating equal powers of (^ + |) we 
arrive at the following expression for the potential 

V(x) = h( \ )x^ + ( )x^ 

^ ' 12 ^2 4 24 48^ ^12 24 144^ 

+(l20-2So'^° + l680^'- 
Thus to lowest order one can think of the Q-oscillator with small values of the pa- 
rameter T as a harmonic oscillator perturbed by a x^ term. A similar expression 
is found, for example, by Taylor expanding the modified Poschl-Teller potential e2, 
which, amouff other applications, has been recently used in the description of hyper- 
nuclei e3e9SE3. The modified Poschl-Teller potential (see also subsec. 32.2) has the 
form 

, , 4 . . 



Its Taylor expansion is 

V{x)pT = Ai-l + - -aV + —a^x"" -...)• (13.10) 
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We remark that this expansion contains the same powers of x as the expansion (13.8). 
Furthermore, the signs of the coefficients of the same powers of x in the two expansions 
are the same for T < 0. 



13.3. WKB-EPs for the q-deformed oscillator 

The potentials obtained above are only rough lowest order estimates. More accu- 
rate methods exist for constructing WKB equivalent potentials (WKB-EPs) giving 
(within the limits of WKB approximation) the same spectrum as the above mentioned 
oscillators. A method by which this can be achieved has been given by Wheeler H and 
is described by Chadan and Chabatier El. Applying this method to the o-deformed 
oscillator (sec. 10) with q being a phase factor one finds the potential 00 
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while for the g-deformed oscillator with q real one has 
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The results of this subsection are corroborated by an independent study of the relation 
between the g-deformedharmonic oscillator and the ordinary anharmonic oscillator 
with x^ anharmonicitie^. 



13.4. WKB-EPs for the Q-deformed oscillator 



Using the same technique one finds that the WKB enuivalent potential for the 
Q-deformed harmonic oscillator (sec. 11) takes the form 153 
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(13.15) 
(13.16) 
(13.17) 



^ ujm{Q + 1) / 

We remark that this WKB-EP contains all even powers of x, in contrast to the 
WKB-EPs for the g-oscillator (eqs (13.11), (13.13)), which contains only the powers 
x^, x^, .... This is in agreement with the lowest order results obtained in subsecs 
13.1 and 13.2. 



14. The quantum algebra suq(2) 

Quantum algebras are generalized versions of the usual Lie algebras, to which they 
reduce when the deformation parameter g is set equal to unity. A simple example of 
a quantum algebra is provided by sug(2) IlltJ, which is generated by the operators 
J+, Jo, J-, satisfying the commutation relations 

[Jo,J±] = ±J±, (14.1) 

[J+,J_] = [2Jo], (14.2) 

with Jg = Jq^ {J+V = We remark that eq. (14.1) is the same as in the case of 
the ordinary su(2) algebra, while eq. (14.2) is different, since in the usual su(2) case 
it reads 

[J+, J_] = 2 Jo. (14.3) 

In the rhs of eq. (14.3) one has the first power of the Jo operator, while in the rhs 
of eq. (14.2) one has the g-operator [2Jo], defined in sec. 2. Because of eqs. (2.2) 
and (2.3) it is clear that if one writes the rhs of eq. (14.2) in expanded form all odd 
powers of Jo will appear: 



Thus sug(2) can be loosely described as a nonlinear generalization of the usual su(2): 
While in usual Lie algebras the commutator of two generators is always producing a 
linear combination of generators, in the case of quantum algebras the commutator of 
two generators can contain higher powers of the generators as well. 

The irreducible representations (irreps) of sUq(2) (which have dimensionality 
2J + 1) are determined by highest weight states with J = 0, |, 1, . . .. The basic states 
|J, M > (with —J < M < J) are connected with highest weight states \J, J > as 
follows 



\J,M >-- 



[J + M]\ 



> J-M 



J, J >, 



(14.6) 



\ [2J]\[J-M]\ 

with J^\J, J >= and < J, J\J, J >= 1. The action of the generators of the algebra 
on these basic vectors is given by 



Jo\J,M>= M\J,M>, 



;i4.7) 



J±\J,M>= ,J[Jt M][J ±M +1]\J,M ±1 > . (14.8) 

These expressions look similar to the ones of the usual su(2) algebra, the only differ- 
ence being that g-numbers appear under the square root instead of ordinary numbers. 

The second order Casimir operator of sug(2) is determined from the condition that 
it should commute with all of the generators of the algebra. The resulting operator 
is 



.2 - J-J+ + [Jo][Jo + 1] = J+J- + [Jo][Jo 

Its eigenvalues in the above mentioned basis are given by 



CI\J,M >= [J][J +1]\J,M >, 



(14.9) 



(14.10) 



while for the usual su(2) the eigenvalues of the Casimir operator are J(J + 1). One 
can easily check that for real q {q = with r real) the eigenvalues [J] [J + 1] produce 
a spectrum increasing more rapidly than J{J + 1) (an expanded spectrum), while for 
q being a phase factor (g = e*'^ with r real) the eigenvalues [J][J + 1] correspond to 
a spectrum increasing less rapidly than J(J + 1) (a compressed spectrum). 

It should be noticed that the generators J+, Jq, J_ of sUg(2) are connected to the 
generators j+, jo, j_ of the usual su(2), which satisfy the commutation relations 



bo, j±] = ±i 



= 2 jo, 



;i4.ii) 



through the relations Sii 



-''o = jo. 



J4 



[jo + j]bo - 1 - j] . 

^(jo+j)(jo-i-jy^' 



[jo + j] [jo 



J] 



(14.12) 

where j is determined by the relation for the second order Casimir operator of su(2) 



C = 3-3+ + Jo(jo + 1) = 3+3- + Jo(jo - 1) = J(J + !)• (14.13) 



15. Realization of sug(2) in terms of g-deformed bosons 



Realizations of Lie algebras in terms of (ordinary) bosons are useful not only as 
a convenient mathematical tool, but also because of their applications in physics 1^. 
In the case of quantum algebras it turns out that boson realizations are possible in 
terms of the g-deformed boson operators already introduced in sec. 10. 

In the case of sUg(2) the generators can be mapped onto g-deformed bosons as 
follows 0'0 ^ 

J+ = a{a2, J- = alai, Jq = -{Ni - N2), (15.1) 

where a|, and Ni are g-deformed boson creation, annihilation and number operators 
as these introduced in sec. 10. One can easily prove that the boson images satisfy 
the commutation relations (14.1) and (14.2). For example, one has 

[J^, J_] = ala^alai - ala,a\a2 = miN^ + 1] - [N, + 1][N2] = [N, - N^] = [2Jo], 

where use of the identity (2.6) has been made. 

In the g-boson picture the normalized highest weight vector is 

|jj>= A^i^lO >, (15.3) 



'[2 J] 

while the general vector \JM > is given by 



( Uj+m ( J\J-m 
\JM >= |0 > . (15.4) 



[J + M]\ ^[J-M]\ 

It should be noticed that it was the search for a boson realization of the sUg(2) 
^gbra that led to the introduction of the g-deformed harmonic oscillator in 1989 

Starting from sUq(2) one can formulate a g-deformed version of angular momentum 
theory. Some references are listed here: 

i) Clebsch-Gordan coeeficients for su,(2) can be found in fflUBiSilii. 

ii) 6-j symbols for sug(2) can be found in @0JlllEi. 

iii) 9-j symbols for sug(2) can be found in 0. 

iv) The g-deformed version of the Wigner-Eckart theorem can be found in 

In addition, it should be noticed that a two-Darameter deformation of su(2), la- 
belled as sup,g(2) has been introduced Clebsch-Gordan coefficients for 
sup,g(2) have been discussed in SSiSH. 

The way in which the algebra sUg(2) can be realized in terms of the Q-deformed 
bosons of sec. 11 is given in subsec. 33.1 (see eqs (33.13)-(33.16)). 



16. The quantum algebra sUq(l,l) 



In this section we shall give a brief account of the algebra suq(l,l) 

In the classical case the so(2,l) generators satisfy the commutation relations E3 

[Ki, K2] = -tKs, [K2, Ks] = iK,, [K3, Ki] = iK2, (16.1) 

which differ from the classical so(3) commutation relations in the sign of the r.h.s. of 
the first commutator. Defining 

K+ = Ki+ K^ = Ki- K3 = K,, (16.2) 

one obtains the su(l,l) commutation relations 

[K,, K±] = ±K±, [K+, K_] = (16.3) 

which differ from the familiar su(2) commutation relations in the sign of the r.h.s. of 
the last commutator. The generators of su(l,l) accept the following boson represen- 
tation H ^ 

= alal, K_ = aia2, = -{a[ai + 0202 + 1), (16.4) 

where a\, ai, al, a2 satisfy usual boson commutation relations. 
The second order Casimir operator of so(2,l) is @ 

C2[so(2, 1)] = -{Kl + Kl - Kl), (16.5) 

while for su(l,l) one has 

C2[su(l, 1)] = [M[K^ - 1] - K+K^ = [K^][K^ + 1] - K^K+. (16.6) 
In the quantum case, the generators of sug(l,l) satisfy the commutation relations 

[K,,K±] = ±K±, [K+,K^] = -[2Ko]. (16.7) 
The generators of sUg(l,l) accept the following boson representation 

K+ = a{al K_ = aiOs, Ko = ^Ni + N2 + 1), (16.8) 

where the bosons aj, ai {i = 1,2) satisfy the usual g-boson commutation relations. 
The second order Casimir operator of sUg(l,l) is 

C2[su,(l, 1)] = [Ko][Ko - 1] - K+K_ = [Ko][Ko + 1] - K_K^. (16.9) 

Its eigenvalues are 

C2[su,(l, l)]|/t/i >= [k][k - l]\KfI >, (16.10) 

where 

l + |ni-n2| l + ni+ns 

K = , n = , (16.11) 



since the basis has the form |fi;yU >= |r2i > |r;,2 >, with 



Ui >= ^^{alp\0 > . (16.12) 



In this basis the possible values of fi are given by 

= K, K + 1, K + 2, . . . , (16.13) 

up to infinity, while k may be any positive real number. The action of the generators 
on this basis is given by 

Ko\Kfi)=fiM, (16.14) 



K±\Kfi) = ^[fi±K\[fiTf^±l]\i^f^±l)- (16.15) 

Clebsch-Gordan and Racah coefficients for sUq(l,l) can be found in HSBBlll. 
Furthemore, a twviiaiajrieter deformed version of sUq(l,l), labelled as sUp^g(l,l), has 



been introduced O'EalHSH 

17. Generalized deformed su(2) algebras 

In the same way that in addition to the g-deformed oscillators one can have 
generalized deformed oscillators, it turns out that generalized deformed su(2) algebras, 
containing sUq(2) as a special case andhaving representation theory similar to that 
of the usual su(2), can be constructed eB'M. It has been proved that it is possible to 
construct an algebra 

[Jo, J±] = ±J±, [J+, J-] = HMJo + 1)) - HMJo - 1)), (17.1) 

where Jo, J+, J- are the generators of the algebra and $(x) is any increasing entire 
function defined for x > —1/4. Since this algebra is characterized by the function $, 
we use for it the symbol su$(2). The appropriate basis \l,m > has the properties 

Jo\L,M >= M\L,M >, (17.2) 



J+|L,M>= J<I)(L(L + 1))-$(M(M + 1))|L,M + 1>, (17.3) 



J_|L,M>= ^<I>(L(L + 1))-$(M(M-1))|L,M-1 >, (17.4) 

where 

L = 0,^,l,|2,|3,..., (17.5) 

and 

M = -L, -L + 1, -L + 2, . . . , L - 2, L - 1, L. (17.6) 
The Casimir operator is 



C = J_J+ + $(Jo(Jo + 1)) = J+J- + $(Jo(Jo - 1)), 



(17.7) 



its eigenvalues indicated by 



C\L,M >=<^{L{L + 1))\L,M > . (17.8) 
The usual su(2) algebra is recovered for 

<^{x{x + 1)) = x{x + 1), (17.9) 
while the quantum algebra sug(2) 

[Jo, J±] = ±J±, [J+, J-] = [2Jo]„ (17.10) 

occurs for 

<^{x{x + 1)) = [x][x + 1]. (17.11) 

The su$(2) algebra occurs in several cases, in which the rhs of the last equation 
in (17.1) is an odd function of Jq eH. It can be seen that other algebraic structures, 
like the quadratic Hahn algebra QH(3) ^ and the finite W algebra Wq ^ can be 
brought into the su$(2) form, the advantage being that the representation theory of 
su<j, (2) is already known. It can also be proved that several physica l syste ms, like the 



isotropic oscillator in a 2-dim curved space with constant curvature I ""11 "H, the Kepler 



system in a 2-dim curved space with constant curvature , and the system of two 
identical particles in two dimensions E^l can also be put into an su$(2) form. More 
details can be found in 0. 

18. Generalized deformed parafermionic oscillators 

It turns out that the generalized deformed su$(2) algebras mentioned in the last 
section are related to generalized deformed parafermionic oscillators, which we will 
therefore describe here. 



It has been proved 11221 that any generalized deformed parafermionic algebra of 



order p can be written as a generalized oscillator (sec. 12) with structure function 

F{x) = x{p + 1 - x){X + fix + vx^ + px^ + ax'^ + .. .), (18.1) 
where A, /i, u, p, a, . . . are real constants satisfying the conditions 

X + px + ux'^ + px^ + ax'^ + . . . > 0, X e {1,2, . . . (18.2) 
Considering an su$ (2) algebra S with structure function 

$(Jo(Jo + 1)) = AJoiJo + 1) + B{Jo{Jo + 1))' + C{Jo{Jo + 1))', (18.3) 
and making the correspondence 

J+ ^ A\ J_ ^A, Jo ^ N, (18.4) 



one finds that the su$(2) algebra is equivalent to a generalized deformed parafermionic 
oscillator of the form 

F{N) = N{p+1- N) 

[-ip\p+ l)C+pB) + (p^C + ip- 1)B)N+ -p+i)C + B)N^ + {p-2)CN'' + CN% 

(18.5) 

if the condition 

A + p{p+l)B + p'^{p+lfC = (18.6) 

holds. The condition of eq. (18.2) is always satisfied for B > and C > 0. 

In the special case of C = one finds that the su$(2) algebra with structure 
function 

$( Jo( Jo + 1)) = AJoiJo + 1) + BiMJo + l)f (18.7) 
is equivalent to a generalized deformed parafermionic oscillator characterized by 

F{N) = BN{p +l-N){-p+{p- 1)N + N^), (18.8) 

if the condition 

A + p{p+l)B = (18.9) 

is satisfied. The condition of eq. (18.2) is satisfied for B > 0. 

Including higher powers of Jo{Jq + 1) in eq. (18.3) results in higher powers of 
in eq. (18.5) and higher powers oi p{p + 1) in eq. (18.6). If, however, one sets B = 
in eq. (18.7), then eq. (18.8) vanishes, indicating that no parafermionic oscillator 
equivalent to the usual su(2) rotator can be constructed. 

It turns out that several other mathematical structures, like the finite W alge- 
bras Wo and W3^'' (see subsec. 34.5) can be put into the generalized deformed 
parafermionic oscillator form. The same is true for several physical systems, such as 
the isotropic os cillator and the Kepler problem in a 2-dim curved space with cnii- 



stant curvature and the Fokas-Lagerstrom t^J, SiQorodinsky-Winternitz 



and Holt t2a potentials. Further details can be found in t3. 
19. The suq(2) rotator model 
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It has been suggested by Raychev, Roussev and Smirnov Eril and independently 



by Iwao I i that rotational spectra of deformed nuclei can be described by the 



g-deformed rotator, which corresponds to the 2nd order Casimir operator of the 
quantum algebra sug(2), already studied in sec. 14. We shall show here that this 
assumption works and discuss the reasons behind this success, as well as the relation 
ti3 between the sUg(2) model and the Variable Moment of Inertia (VMI) model (see 
sec. 20). 

The g-deformed rotator corresponds to the Hamiltonian 



(19.1) 



where / is the moment of inertia and Eo is the bandhead energy (for ground state 
bands Eq = 0). For q real, i.e. with q = with r real, the energy levels of the 
g-rotator are 



EiJ) ^ + 1] + ^0 = ^^''^Z^ + + E„. (19.2) 

For q being a phase, i.e. q = e*"^ with r real, one obtains 



Raychev et al. 1123 have found that good fits of rotational spectra of even-even 
rare earths and actinides are obtained with eq. (19.3). It is easy to check that eq. 

(19.2) fails in describing such spectra. In order to understand this difference, it is 
useful to make Taylor expansions of the quantities in the numerator of eq. (19.2) (eq. 

(19.3) ) and collect together the terms containing the same powers of J(J + 1) (all 
other terms cancel out), finally summing up the coefficients of each power. In the 
first case the final result is 

J( J + 1))' + Y^I]^'v2(.r)(AJ + l))-" + . . . (19.4) 

where (r) are the modified spherical Bessel functions of the first kind 0. 

In the second case (eq. (19.3)) following the same procedure one obtains 

EiJ) = Eo + ^ (-^Oo(r) J(J + 1) - rnir)iJiJ + 1))^ 

+ lr'Mr){J{J + 1))' - lr'MT){J{J + 1))' + ^r%{r){j{j + 1))^ -...), (19.5) 

where jn{T) are the spherical Bessel functions of the first kind IHl. 
Both results are of the form 

E{J) = AJ{J + 1) + B{J{J + 1))2 + C( J( J + 1))3 + D{J{J + 1))^ + . . . . (19.6) 

Empirically it is known that nuclear rotational spectra do show such a behaviour, the 
coefficients A, B, C , D, . . . having alternating signs (starting with A positive) and 
magnitudes dropping by about 3 orders of magnitude each time one moves to the 
next higher power of J( J + 1) Eil'Eil. 

It is interesting to check if the empirical characteristics of the coefficients A, B, 
C, D are present in the case of the expansions of eqs. (19.2), (19.3), especially for 
small values of r. (Since we deal with rotational spectra, which are in first order 



approximation described by the usual algebra su(2), we expect r to be relatively 
small, i.e. the deviation of sUq(2) from su(2) to be small. This is in agreement to the 
findings of t23, where r is found to be around 0.03.) 

One can easily see that in eq. (19.2) it is impossible to get alternating signs, while 
in eq. (19.3) the condition of alternating signs is readily fulfilled. This fact has as 
a result that the energy levels given by eq. (19.2) increase more rapidly than the 
levels given by the J( J + 1) rule, while the levels given by eq. (19.3) increase less 
rapidly than J {J + 1). In order to check the order of magnitude of the coefficients 
for small values of r, it is useful to expand the spherical Bessel functions appearing 
in eq. (19.3) and keep only the lowest order term in each expansion. The result is 

E{j) =E,+ + 1) - + i)r + + i)r 

+ i)r + -^{JiJ + 1))' -•••)• (19.7) 
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We remark that each term contains a factor more than the previous one. For r in 
the area of 0.03, is of the order of 10~^, as it should. We conclude therefore that 
eq. (19.6) is suitable for fitting rotational spectra, since its coefficients have the same 
characteristics as the empirical coefficients of eq. (19.6). 

Extended comparisons of the sUg(2) predictions to experimental data for ground 
state bands of rare earth and actinide nuclei can be found in Es^Slil^lili. More 
recently, the sUg(2) formalism has been used for the description of /3- and 7-bands of 
deformed rare earths and actinides, with satisfactory results clB. 

It is necessary for E{J) to be an increasing function of J. In order to guarantee 
this in eq. (19.3) one must have 

r(J+l)<^. (19.8) 



In the case of r = 0.036 (as in ^sa^j |iio|-j^ finds J < 42, this limiting value being 
larger than the highest observed J in ground state bands in the actinide region 1113. 
Similarly, for r = 0.046 (as in ^^^Hf in 11121), one finds J < 32, this limiting value being 



again higher than the highest observed J in ground state bands in the rare earth 
region ti3. 

20. Comparison of the suy(2) model to other models 



20.1. The Variable Moment of Inertia (VMI) model 

In lowest order approximation rotational nuclear spectra can be described by the 
formula 

EiJ) = (20.1) 



where B is the moment of inertia of the nucleus, which is assumed to be constant. 
However, in order to get closer agreement to experimental data, one finds that he has 
to include higher order terms in this formula, as shown in eq. (19.6). 

Another way to improve agreement with experiment is to let the moment of inertia 
O to vary as a function of the a ngu lar momentum J. One thus obtains the Variable 
Moment of Inertia (VMI) model In this model the levels of the ground state band 
are given by 

E{J) = + \cmJ) - eo)^ (20.2) 

where 0(J) is the moment of inertia of the nucleus at the state with angular momen- 
tum J, while C and Bq are the two free parameters of the model, fitted to the data. 
The parameter 0o corresponds to the ground state moment of inertia, while instead 
of the parameter C it has been found meaningful to use the parameter combination 

" = del- 

which is related to the softness of the nucleus. The moment of inertia at given J is 
determined through the variational condition 

which is equivalent to the cubic equation 

^3 o.n2o ^(•^+1) 



Q{JY - QiJyeo - ^-1^ = 0- (20.5) 



This equation has only one real root, which can be written as 

e(j) 



j(j + i) ei l {J{J + i)y egj(j + i) 
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J(J+1) , eg (J(J+1))2 , QpiJ+l) , Go 

Expanding the roots in this expression one obtains 

e( J) = 0o(l + aJ( J + 1) - 2a\j{J + 1))2 

+7a\j{J + l)f - 30a\JiJ + 1))^ + ...). (20.7) 
Using eq. (20.7) in eq. (20.2) one obtains the following expansion for the energy 

E{J) = ^^{J{J + l)-lcr{J{J+l)r 



W{J{J + 1))=^ - 3a3( J( J + l)Y + ...). (20.8) 

Empirically it is known 111! that for rotational nuclei the softness parameter a is of the 
order of 10~^. Therefore the expansion of eq. (20.8) has the same characteristics as 
the expansion of eq. (19.6) (alternating signs, successive coefficients falling by about 
3 orders of magnitude). 



20.2. Comparison of the suq(2) model to the VMI and related models 

We now turn to the comparison of the expansion of eq. (19.5) to the Variable 
Moment of Inertia (VMI) model, discussed in the previous subsection. Comparing 
eqs (19.5) and (20.8) we see that both expansions have the same form. The moment 
of inertia parameter / of (19.5) corresponds to the ground state moment of inertia 
00 of (20.8). The small parameter of the expansion is in the first case, while it is 
the softness parameter l/(2C0o) in the second. However, the numerical coefficients 
in front of each power of J( J + 1) are not the same. 

In a comparison is made between the parameters obtained by fitting the same 
spectra by the sug(2) and VMI formulae. The agreement between 1/(2/) and l/(20o) 
is very good, as it is the agreement between and a as well. Therefore the known 
Ell smooth variation of 0o and a with the ratio = E^A) L£li2) is expected to hold 
for the parameters / and as well. This is indeed seen in 

The difference between the expansions of eqs (19.5) (or (19.7)) and (20.8) is also 
demonstrated by forming the dimensionless ratios AC/{AB'^) and A^D/{24:B^) in eq. 
(19.6). In the case of eq. (20.8) bothgumitities are equal to 1, as expected, since it 
is known that the VMI is equivalent tl^'EHj to the Harris expansion ll2iJ, in which both 
quantities are known to be equal to 1. In the case of eq. (19.7) the corresponding 
values are AC/{AB'^) = 1/10 and A^D/{2AB^) = 1/280. According to the empirical 
values of these ratios given in llHllli, the ratios given from eq. (19.7) are better than 
the ratios given by eq. (20.8), especially the second one. 



20.3. The hybrid model 



The hybrid model of nuclear collective motion 



123,124,125 



'E^ has been introduced in 

order to provide a link between the two successful ways of describing lo w-lying nuclear 
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tHa and the 
The hybrid 



excitations: the extended form of the Bohr-Mottelson model (BMM) 
Interacting Boson Model (IBM) tH3 (see sees 27, 29 for more details), 
model combines the advantages of both models, i.e. the geometrical significance of the 
collective coordinates inherent in the extended BMM, and the use of group theoretical 



concepts characterizing IBM. In the framework of the rotational limi^f t 

r 



le hybrid 

model, associated with the chain u(6)Dsu(3), Partensky and Quesne ll22Hli!ll, starting 
from the fact that in the geometrical description the square of the deformation is 
proportional to the moment of inertia of the ground state band, proved that the 



energy levels of the ground state band are given by 



with A being a free parameter and B given by 

5 = 8A^2 + 22A^ - 15, (20.10) 

where N is the sum of the number of valence proton pairs (or proton-hole pairs, when 
more than half of the proton valence shell is filled) A^^r and the number of the valence 
neutron pairs (or neutron-hole pairs, when more than half of the neutron valence shell 
is filled) Njj. We remark that in the framework of this model the moment of inertia 

e«A.).^(^±i)±^ (20.11) 

depends on both the angular momentum J and the valence pair number A^. 
It is instructive to expand E{J) in powers of J(J+1) t23 

~{J{J + 1))' + + 1))' - + 1))' + •••)• (20.12) 

Comparing the present expansion to the one of eq. (19.7) for the sUg(2) model one 
has 

'3 / 3 




w o , (20.13) 

B V 8iV^ + 22iV - 15 ^ ^ 

obtaining thus a connection between the r parameter of the suq(2) model and a 
microscopic quantity, the valence nucleon pair number. From this equation it is clear 
that T decreases for increasing N. Thus the minimum values of r are expected near 
the midshell regions, where the best rotators are known to be located. This is a 
reasonable result, since r describes the deviations from the su(2) (rigid rotator) limit. 
It is clear that the minimum deviations should occur in the case of the best rigid 
rotators. 

From eq. (20.13) one can obtain for each nucleus the value of r from the number 
of valence nucleon pairs present in it. These predictions for r have been compared 
in llH to the values of r found empirically by fitting the corresponding spectra, with 
good results in both the rare earth and the actinide regions. In addition eq. (20.13) 
indicates that in a given shell nuclei characterized by the same valence nucleon pair 
number will correspond to the same value of r. Such multiplets have been studied 
in the framework of the hybrid model by EH. 

It is worth noticing that fits of 71-bands in the rare earth region (Er, Yb isotopes) 
give T parameters very similar to the ones coming from fitting the corresponding 
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ground state bands, in addition exhibiting the same as the one mentioned above 
behaviour of r as a function of A^. 

Taking further advantage of the above connection between the sUq(2) model and 
the hybrid model, the parameter r has been connected to the nuclear deformation 
parameter /3, as well as to the electromagnetic transition probabilities B(E2:25'" —>■ Of). 
Since both f3 and B(E2:2^ Of) are known to increase with increasing collectivity, 
r is expected to decrease with increasing f3 or increasing B(E2:2^ Of). This 
expectation is corroborated by the results reported in 
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20.4- Other models 

It should be noticed that an empirical formula very similar to that of eq. (19.3) 
has been recently proposed by Amal'skii 1123 on completely different physical grounds. 
The formula reads 

2 /TtJ 



E{J)=Asm'^—j, (20.14) 

where A and B are free parameters. This formula should be compared to eq. (19.3), 
with which it is almost identical. 

A different formula, also giving very good fits of rotational spectra, has been 



introduced by Celeghini, Giachetti, Sorace and Tarlini IHj, based on the g-Poincare 
rotator. 



21. Electromagnetic transitions in the su^(2) model 

We have already seen that the sug(2) formalism provides an alternative to the 
VMI model, the deformation parameter q being connected to the softness parameter 
of the VMI model. 

The stretching effect present in rotational energy levels, which can equally well be 
described in terms of the VMI model and the sUg(2) symmetry, should also manifest 
itself in the B(E2) transition probabilities among these levels. If deviations from 
the su(2) symmetry are observed in the energy levels of a band, relevant deviations 
should also appear in the B(E2) transitions connecting them. In the case of the 
VMI model no way has been found for making predictions for the B(E2) transition 
probabilities connecting the levels of a collective band. The sug(2) symmetry naturally 
provides such a link. Before studying the suq(2) case, though, it is useful to recall the 
predictions of other models on this matter. 



21.1. The collective model of Bohr and Mottelson 



In rotational bands one has 11^ 

B{E2 : J + 2 J) = -^Ql\Cifo':^'j^\^ (21.1) 



i.e. the B(E2) transition probability depends on the relevant Clebsch-Gordan coef- 
ficient of su(2), while Qq is the intrinsic electric quadrupole moment and K is the 
projection of the angular momentum J on the symmetry axis of the nucleusin, the 
body-fixed frame. For K = bands, as the ground state bands, one then has 
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B(£2:J + 2^J) = AqJ? (£±1)(£±^, (21.2) 

^ ^ 167r^°2(2J + 3)(2J + 5) ^ ^ 

It is clear that the B(E2) values should saturate with increasing J. 

21.2. The Interacting Boson Model (IBM) 

It is also instructive to mention what happens in the case of the Interacting Boson 
Model (IBM) llHa, the successful algebraic model of nuclear structure with which we 
are going to deal in more detail later (see sections 27, 29). In the case of the su(3) 
limit of the IBM, which is the limit applicable to deformed nuclei, the corresponding 
expression is El^ 

B(E2 :J^2^J) = Aggg (J^l)(J^2) (2N- J)(2N + J + 3) 

^ ' 167r^°2(2J + 3)(2J + 5) (2A^ + 3/2)2 ' v ; 

where N is the total number of bosons. Instead of saturation one then gets a decrease 
of the B(E2) values at high J, which finally reach zero at J = 2N. This is a well 
known disadvantage of the simplest version of the model (IBM-1) due to the small 
number of collective bosons {s {J = 0) and d [J = 2)) taken into account. It can 



be corrected by the inclusion of highe 
approximately restore saturation (see E 



sons {g {J = A), i {J = Q), etc), which 
for full list of references). 



Another way to avoid the problem of decreasing B(E2)s in the su(3) limit of IBM 
at high J is the recently proposed 1123 transition from the compact su(3) algebra to 
the noncompact sl(3,R) algebra. The angular momentum at which this transition 
takes place is fitted to experiment. In this way an increase of the B(E2) values at 
high J is predicted, which agrees well with the experimental data for ^^^U. 

21.3. The suq(2) model 

In order to derive a formula similar to (21.2) in the sug(2) case, one needs to 
develop an sug(2) angular momentum theory. As mentioned in sec. 15, this has 
already been achieved. It turns out that an equation similar to (21.1) holds in the 
q-deformed case, the only difference being that the Clebsch-Gordan coefficient of the 
sUg(2) algebra must be used instead. These coefficients have the form E3 



^J+2,2,J _ 2K 
q^K,0,-K — H 



\^\A][J - K + 2][J - K + l][J + K + l][J + K + 2] 
\\ [2] [2 J + 2] [2 J + 3] [2 J + 4] [2 J + 5] ' ^ ' ^ 



For K = bands one then has 



' ' 167r^°[2][2J + 2][2J + 3][2J + 4][2J + 5] ^ ' 

For q = e*"^ this equation takes the form 

,sin(3r) sin(4r) 



B,{E2;J + 2^J) = —Q',- 



IGtt sin(2r) sin(r) 

(sin(r(J + l)))^(sin(r(J + 2)))^ 

sin(r(2J + 2))sin(r(2J + 3))sin(r(2J + 4))sin(r(2J + 5))' ^ ' ^ 

It is useful to get an idea of the behaviour of this expression as a function of J, 
especially for the small values of r found appropriate for the description of ground 
state spectra. Expanding all functions and keeping corrections of the leading order 
in r only, one has 



5 (J+l)(J + 2) 



r 



2 



B,(E2, J + 2 ^ J) = + -(U' + 22J + 12)), (21,7) 

We see that the extra factor, which depends on r^, contributes an extra increase 
with J, while the usual su(2) expression reaches saturation at high J and IBM even 
predicts a decrease. 

21.4- Comparison to experiment 

Is there any experimental evidence for such an increase? In order to answer this 
question one should discover cases in which the data will be consistent with the sug(2) 
expression but inconsistent with the classical su(2) expression. (The opposite cannot 
happen, since the classical expression is obtained from the quantum expression for 
the special parameter value r = 0.) Since error bars of B(E2) values are usually large, 
in most cases both symmetries are consistent with the data. One should expect the 
differences to show up more clearly in two cases: 

i) In rare earth nuclei not very much deformed (i.e. with an = E{A) / E{2) ratio 
around 3.0). These should be deformed enough so that the sug(2) symmetry will be 
able to describe them, having, however, at the same time values of r not very small. 
Since in several of these nuclei backbending (or upbending) occurs at J = 14 or 16, 
one can expect only 5 or 6 experimental points to compare the theoretical predictions 
with. 

ii) In the actinide region no backbending occurs up to around J = 30, so that 
this is a better test ground for the two symmetries. However, most nuclei in this 
region are well deformed, so that small values of r should be expected, making the 
distinction between the two theoretical predictions difficult. 



A few characteristic examples (4 rare earths and an actinide) are given in 11121. In 



all cases it is clear that the suq(2) curve follows the experimental points, while the 



su(2) curve has a different shape which cannot be forced to go through all the error 
bars. Several comments are now in place: 

i) For a given nucleus the value of the parameter r obtained from fitting the 
B(E2) values among the levels of the ground state band should be oqual to the value 
obtained from fitting the energy levels of the ground state band. In tf3 it is clear that 
both values are similar, although in most cases the value obtained from the B(E2)s 
is smaller than the value obtained from the spectra. It should be taken into account, 
however, that in most cases the number n' of levels fitted is different (larger) than 
the number n of the B(E2) values fitted. In the single case (^®^W) in which n = n', 
the two r values are almost identical, as they should. 

ii) One can certainly try different fitting procedures. Using the value of r obtained 
from the B(E2) values for fitting the spectrum one gets a reasonably good description 
of it, although the squeezing of the spectrum is not as much as it should have been 
(with the exception of ^^^W). Using the value of r obtained from the spectrum for 
fitting the B(E2) values one obtains an increase more rapid than the one shown by 
the data (again with the exception of ^^"^W). One can also try to make an overall fit 
of spectra and B(E2)s using a common value of r. Then both the squeezing of the 
spectrum and the rise of the B(E2)s can be accounted for reasonably well although 
not exactly. One should notice, however, that the experimental uncertainties of the 
B(E2)s are much higher than the uncertainties of the energy levels. 

iii) Concerning energy levels, the rigid rotator model and the su(3) limit of the 
IBM predict a J(J+ 1) increase, while the sUg(2) model and the VMI model predict 
squeezing, which is seen experimentally. 

iv) Concerning the B(E2) values, the VMI makes no prediction, the rigid rotator 
predicts saturation at high J, the su(3) limit of the IBM predicts decrease, while the 
sug(2) model predicts an increase. The evidence presented in tlj supports the sug(2) 
prediction, but clearly much more work, both experimental and analytical, is needed 
before final conclusions can be drawn. The modified su(3) limit of IBM described in 

also supports the increase of the B(E2) values at high J. Increasing BE(2) values 
are also predicted in the framework of the Fermion Dynamical Symmetry Model 
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There is also empirical evidence for increasing B(E2) values in the recent systematics 
by Zamfir and Casten llil. 

v) It is clear that much further work is needed as far as comparisons of the sUg(2) 
predictions to experimental BE(2) values are concerned for safe conclusions to be 
reached. 

vi) Since the quadrupole operator is not a member of the symmetry algebra sug(2) 
under consideration, it is clear that the B(E2) values studied here do not contain any 
dynamical deformation effects, but only the kinematical ones (through the use of the 
g-deformed Clebsch-Gordan coefficients). A more complete approach to the problem 
will be the construction of a larger algebra, of which the quadrupole operator will be 
a member and it will also be an irreducible tensor under sUq(2) or sOq(3). Work in 
this direction is in progress. 



22. Superdeformed bands 



One of the most impressive experimental disc overies in nuclear physics during the 
last decade was that of superdeformation S (see for relevant reviews). The 

energy levels of superdeformed bands follow the J{J+ 1) rule much more closely than 
the usual rotational bands. Levels with J larger than 60 have been observed. A 
compilation of superdeformed bands has been given in llH. The best examples have 
been found in the A^150 mass region, while additional examples have been found 
in the A^80, A^130 and A^190 regions. It is understood that the superdeformed 
bands in the A^150 region correspond to elongated ellipsoidal shapes with an axis 
ratio close to 2:1, while in the A^ISO and A^190 regions the ratios 3:2 and 1.65:1 
respectively appear closer to reality. 

Since the sUg(2) model has been found suitable for describing normal deformed 
bands, it is plausible that it will also be successful in describing superdeformed bands 



as well. A test has been performed in Ilia. The sUg(2) model has been found to give 



good results in all mass regions, the deformation parameter r being smaller than in 
the case of normal deformed bands, thus indicating smaller deviations from the sug(2) 
symmetry. In particular, r has been found to obtain values about 0.01 in the A;^130 
and A:^190 regions, while it obtains even smaller values, around 0.004, in the A;^150 
region, which contains the best examples of superdeformed bands observed so far. 
These results should be compared to the values of r around 0.03 obtained in the case 
of normal deformations. 

Concerning the corresponding B(E2) values, the experimental information is still 
quite poor for allowing a meaningful comparison of the sUg(2) predictions to experi- 
ment. 



23. The physical content of the sUq(2) model 

From the above it is clear that the sug(2) model offers a way of describing nuclear 
stretching, i.e. the departure of deformed nuclei from the su(2) symmetry of the rigid 
rotator, similar to the one of the VMI model. The parameter r describes this depar- 
ture quantitatively, vanishing in the rigid rotator limit. Therefore the deformation 
parameter r should not be confused with nuclear deformation; it is in fact related to 
nuclear softness, as already discussed in sec. 20. 

On the other hand, the increase of the moment of inertia with increasing J in the 
framework of the VMI model means that collectivity gets increased ll^. The sug(2) 
model is an alternative way for describing this increase in collectivity. But increased 
collectivity implies increased B(E2) transitions. Therefore it is not surprising that 
the sUg{2) model predicts B(E2) values increasing with J. 

Given the sUg(2) generators J+, J_, Jq, it is instructive to define as usual the 
operators J^, Jy, Jz by 

J+ = Jx ~\~ iJy: J- = Jx ~ iJy, Jq = Jz- (^3.1) 



The suq(2) commutation relations can then be rewritten in the form 

\Jx^ Jyl 2 t^^'^]' I'^y ^Jxi \Jzi Jx\ iJyi (23.2) 

which is a generahzation of the so (3) commutation relations, obtained in the limit 
g ^ 1. We remark that while in the classical so (3) case the three commutation 
relations have exactly the same form, in the quantum case the first commutation 
relation differs (in the right hand side) from the other two, thus indicating that in 
the framework of the problem under study the z-direction is not any more equivalent 
to the X- and y- directions. This is of course a phenomenological way to describe the 
softness of deformed nuclei by adding the appropriate perturbations to the pum^u(2) 
Hamiltonian and has nothing to do with the isotropy of space, as implied in 

It is worth remarking at this point that the Casimir operator of sUg(2) is also 
invariant under the usual su(2) li^. Therefore the quantum number J characterizing 
the irreps of suq(2), and as a result the nuclear levels, is exactly the same as the 
quantum number J used in the case of the usual su(2). Therefore there is no reason for 
a "total reformulation of quantum mechanics", as implied int^Hl, the sug(2) generators 
being connected to their su(2) counterparts by the g-deforming functionals of sec. 14. 
In other words, one continues to believe in usual angular momentum theory and 
usual quantum mechanics. All what is done in the framework of the sUg(2) model 
is to add to the usual su(2) Hamiltonian several perturbations which have a special 
form making them suitable to be summed up, including the original su(2) term, into 
the form of the sUq(2) Hamiltonian. 



24. The Upq(2) rotator model 

An extension of the sug(2) model is the Up^g(2) model lllillllillli, which is based on 
a two-parameter deformed algebra (see sec. 15 for a list of references). For p = q 
(using the definition of {p, g)-numbers of eq. (2.11)) this model reduces to the sU q(2) 
one. This model has been successfully applied to superdeformed nuclear bands 
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When Taylor expanded, it becomes clear that the eigenvalues of the Hamiltonian of 
this model (which is the second order Casimir operator of Up ,j(2)) contain terms of 
the form J( J( J + 1))", in addition to the ( J( J + 1))" ones. It is therefore closer to 
the modification of sug(2) which will be discussed in subsec. 26.3. 



25. Generalized deformed su(2) models 



Another formula giving ver\ 
duced by Holmberg and Lipas I 
levels are given by 

E{J) 



jood results for rotational spectra has been intro- 



and rediscovered in 11^. In this case the energy 



1 + 6J(J + 1) - 1 



(25.1) 



Taylor expansion of the square root immediately shows that the present formula is a 



special case of eq. ( 1 This formula can be derived from the collective model of 

Bohr and Mottelson t2a. 

It has been argued in Ii23 that the Hamiltonian of eq. (25.1) gives better agreement 
to rotational nuclear spectra than the one coming from the suq(2) symmetry. Using 
the techniques described in detail in sec. 17 one can construct a generalized deformed 
algebra su$(2), characterized by a function $(J(J + 1), giving the spectrum of eq. 
(25.1) exactly. In this particular case the algebra is characterized by the structure 
function 

1))=' 



1 + 6J(J + 1)-1 



(25.2) 



It is of interest to check if this choice of structure function also improves the agree- 
ment between theory and experiment in the case of the electromagnetic transition 
probabilities connecting these energy levels. In order to study this problem, one has 
to construct the relevant generalized Clebsch-Gordan coefficients. This problem is 
still open. 



26. Quantum algebraic description of vibrational and transitional nuclear 
spectra 

We have already seen that the sug(2) model describes successfully deformed and 
superdeformed bands. It is not surprising that the applicability of the sug(2) for- 
malism is limited to the rotational region (where the ratio R4 = E{4)/ E{2) obtains 
values between 3.0 and 3.33), since it is based on a deformation of the rotation algebra. 
For describing nuclear spectra in the vibrational (2.0 < -R4 < 2.4) and transitional 
(2.4 < i?4 < 3.0) regions it is clear that an extension of the model is needed. In order 
to be guided towards such an extension, we briefly review the existing experience of 
other successful models. 



26.1. The Interacting Boson Model 

In the rotational (su(3)) limit of the Interacting Boson Model (IBM) (see sees 
27, 29 for more details) the spectrum is described iiv a J{J + 1) expression, while 
in the vibrational (u(5)) and transitional (0(6)) limits expressions of the form 
J( J + c) with c > 1 appear. In the u(5) limit, in particular, the energy levels are 
given by 

E{N, rid, V, riA, J, Mj) = Eq + en^ + aridirid + 4) + f32v{v + 3) + j2J{J + 1), (26.1) 

where is the total number of bosons, is the number of d-bosons, v is the seniority, 
ua is the "missing" quantum number in the reduction from o(5) to o(3), Mj is the 
third component of the angular momentum J, while i?o, e, «, /?, 7 are free parameters. 
The ground state band, in particular, is characterized by quantum numbers = 0, 
1, 2, . . . , V = rid, nj\ = 0, J = 2nd, so that the energy expression for it reads 

E{J) = Eo + ^J + jJ{J + 8) + ^JiJ + 6) + 2^J{J + 1). (26.2) 



In the 0(6) limit the energy is given by 



E{N, a, r, i/a, J, Mj) = Eq + /32r(r + 3) + 72J( J + 1) + r]2cr{a + 4), (26.3) 

where a is the quantum number characterizing the irreducible representations (irreps) 
of 0(6), T is the quantum number characterizing the irreps of o(5), z/a is the missing 
quantum number in the reduction from o(5) to o(3), while Eq, jS, 7, rj are free pa- 
rameters. The ground state band is characterized by the quantum numbers a = N, 
r = 0, 1, 2, . . . , z/A = 0, J = 2r, so that the relevant energy expression takes the 
form 

E{J) = Eo + ^J{J + 6) + 72J(J + 1) + ri2N{N + 4). (26.4) 

The message from eqs (26.2) and (26.4) is that nuclear anharmonicities are de- 
scribed by expressions in which J and appear with different coefficients, and not 
with the same coefficient as in J{J + 1). The earliest introduction of this idea is in 
fact the Ejiri formula 

E{J) = kJ{J +l) + aJ, (26.5) 
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which has been subsequently justified microscopically in 

26.2. Generalized VMI 

The two-parameter VMI model is known to continue giving good fits in the tran- 
sitional and even in the vibrational region. In these regions, however, the accuracy of 
the model is substantially improved by adding a third ^pf^ p.T p eter. which essentially 
allows for treating J and with a different coefficient 
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The usual VMI model has been briefly reviewed in subsec. 20.1. One of the 
(essentially equivalent) three-parameter extensions of the model, which give improved 
fits of vibrational and transitional spectra, is the generalized VMI (GVMI) model 
t^'E^ in which the energy levels are described by 

E{J) = ^^^^^j^'^^ + IhHJ) - $o)^ (26.6) 
which can be easily rewritten in the form 

where x' = — 2. It is clear that for a; = 1/3 the GVMI reduces to the usual 
VMI, while for transitional and vibrational nuclei x obtains lower values, so that x' 
becomes greater than 1. The variational condition determining the moment of inertia 
still has the form of eq. (20.4), while the expansion of the energy turns out to be 

EiJ) = ^y{j+x')-^{j{j+x')r 



+ ((t')'( J( J + x')f - 3{a'f{J{J + x')Y + ...), (26.8) 



where 

We remark that an expansion in terms of J{J + x') is obtained, as compared to an 
expansion in terms of J( J+1) in the case of the usual VMI. The physical content of the 
parameters is clear: the centrifugal stretching effect is accounted for by the softness 
parameter a', as in the case of the usual VMI, while anharmonicities, important in 
the vibrational region, are introduced by x' > 1. Since centrifugal stretching and 
anharmonicities are two effects of different origins, it is reasonable to describe them 
by two different parameters. 

26.3. Modification of the sUq(2) model 

The evidence coming from the IBM and the generalized VMI model described 
above, suggests a model in which the spectrum is given by 

E{J) = + (26.10) 

which contains 3 parameters: the moment of inertia /, the deformation parameter 
q and the new parameter c, which is expected to be 1 in the rotational limit and 
larger than 1 in the vibrational and transitional regions. This energy expression can 
be expanded as 

E{J) = ^ (-^Oo(r) J(J + c) - rn{r){J{J + c)f 

+\r'j2{r){J{J + c)f - \T'n{r){J{J + c)Y + l-^T'u{r){J{J + c)f -...), (26.11) 

which is similar to eq. (19.5) with J(J + 1) replaced by J(J + c). 

It is expected that the deformation parameter r, which plays the role of the small 
parameter in the expansion, as the softness parameter does in the case of the VMI, 
will describe the centrifugal stretching effect, while the parameter c will correspond 
to the anharmonicity effects. These expectations are corroborated from fits of the 



experimental data reported in 11521. The following comments can be made: 

i) The anharmonicity parameter c is clearly decreasing with increasing i?4, i.e. 
with increasing collectivity. It obtains high values (8-18) in the vibrational region, 
while in the rotational region it stays close to 1. (It should be noted that by fixing 
c = 1 in the rotational region the fits are only very slightly changed, as expected.) In 
the transitional region its values are close to 3. 

ii) The deformation parameter r, which corresponds to the centrifugal stretching, 
is known from the sUq(2) model to obtain values close to 0.3-0.4 in the rotational 
region, a fact also seen here. The same range of values appears in the vibrational 



region as well, while in the transitional region r reaches values as high as 0.6. It is not 
unreasonable for this parameter, which is connected to the softness of the nucleus, to 
obtain its highest values in the region of 7-soft nuclei. 

iii) It is worth remarking that eq. (26.10) coincides for q = I and c=integer=N 
with the eigenvalues of the Casimir operator of the algebra o(N+2) in completely 



symmetric states PT l ^^i . In the rotational region the fits gave N=l, which corresponds 
to 0(3), as expected, while in the transitional region the fits gave approximately N=3, 
which corresponds to o(5), which is a subalgebra contained in both the u(5) and 0(6) 
limits of the IBM. 

iv) It is also worth remarking that a special case of the expression of eq. (26.10) 
occurs in the q-deformed version of the 0(6) limit of the Interacting Boson Model, 
which will be reported below (see eq. (29.7)). 

v) The sUg(2) symmetry is known to make specific predictions for the deviation of 
the behaviour of the B(E2) values from the rigid rotator model (subsec. 21.3). It will 
be interesting to connect the spectrum of eq. (26.10) to some deformed symmetry, at 
least for special values of c, and examine the implications of such a symmetry for the 
B(E2) values. Such a study in the framework of the q-deformed version of the 0(6) 
limit of IBM, mentioned in iv), is also of interest. 

vi) It is worth noticing that an expansion in terms of J(J + c) can also be obtained 
from a generalized oscillator (sec. 12) with a structure function 

FiJ) = [JiJ + c)]q, (26.12) 

where [x]q stands for the Q-numbers introduced in sec. 6 and Q = e^, with T real. 
This is similar to an oscillator successfully used for the description of vibrational 
spectra of diatomic molecules (see sec. 35). It can also be considered as a deformation 
of the oscillator corresponding to the Morse potential (see sec. 35). 

We have therefore introduced an extension of the sUq(2) model of rotational nu- 
clear spectra, which is applicable in the vibrational and transitional regions as well. 
This extension is in agreement with the Interacting Boson Model and the Generalized 
Variable Moment of Inertia model. In addition to the overall scale parameter, the 
model contains two parameters, one related to the centrifugal stretching and another 
related to nuclear anharmonicities. In the rotational region the model coincides with 
the usual suq(2) model, while in the transitional region an approximate o(5) symme- 
try is seen. These results give additional motivation in pursueing the construction of 
a deformed version of the Interacting Boson Model. This problem will be discussed 
in the next sections. 



27. A toy Interacting Boson Model with suq(3) symmetry 

The Interacting Boson Model (IBM) EhI is a very popular algebraic model of 
nuclear structure. In the simplest version of IBM low lying collective nuclear spectra 
are described in terms of s (J = 0) and d {J = 2) bosons, which are supposed 
to be correlated fermion pairs. The symmetry of the simplest version of the model 



is u(6), which contains u(5) (vibrational), su(3) (rotational) and 0(6) (7-unstable) 
chains of subalgebras (see also sec. 29). A simplified version of the model, having 
the su(3) symmetry with su(2) and so(3) chains of subalgebras also exists li^. It 
can be considered as a toy model for two-dimensional nuclei, but it is very useful in 
demonstrating the basic techniques used in the full IBM. 

In the present section we will construct the g-deformed version of this toy model. 
Since this project requires the construction of a realization of sUg(3) in terms of q- 
deformed bosons, we will also use this opportunity to study sUg(3) in some detail. 



27.1. The suq(3) algebra 



In the classical version of the toy IBM li^ one introduces bosons with angular 



momentum m = 0, ±2, represented by the creation (annihilation) operators aj, a!^, 
aL (ao, a_|_, a J). They satisfy usual boson commutation relations 

[ai,a\] = 6ij, [ai,aj] = [al,a]] = 0. (27.1) 

The 9 bilinear operators 

Aij = a\aj (27.2) 
satisfy then the commutation relations 

[Aij, Afci] = 6jkKi - ^iAkj, (27.3) 

which are the standard u(3) commutation relations. The total number of bosons 

= SjAjj = ajao + a\a^ + aLa_ (27.4) 

is kept constant. Since we are dealing with a system of bosons, only the totally 
symmetric irreps {A^, 0, 0} of u(3) occur. 

In the quantum case one has the Uq(3) commutation relations given in Table 2 



where A^j are the generators of Ug(3) and the g-commutator is defined as 

[A,B]q = AB -qBA. (27.5) 

In order to obtain a realization of Uq(3) in terms of the g-bosons described in sec. 
10, one starts with 11^ 

A12 = a{a2, A21 = alai, A23 = ala^, A32 = 4«2- (27.6) 

One can easily verify that the Ug(3) commutation relations involving these generators 
are satisfied. For example, one has 

[A12, A21] = [N, - N2], [A23, A32] = [N2 - N,], (27.7) 

using the identity of eq. (2.6) and the identifications 

Ni = Au, N2 = A22, Ns = As3. (27.8) 



Table 2: Ug(3) commutation relations given in the form [A, B]a = C. A is given 
in the first column, B in the first row. C is given at the intersection of the row 
containing A with the column containing B. a, when different from 1, follows C, 
enclosed in parentheses. 
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One can now determine the boson realizations of A13 and A31 from other commutation 
relations, as follows 



Ai3 = [Ai2,A23]g = 
A31 = [A32,A2l]g-i 

Using eq. (2.6) once more one can verify that the relation 

[Ai3,A3i] = [iVi-iV3] 



(27.9) 
(27.10) 

(27.11) 



is fulfilled by the boson images of (27.9), (27.10). It is by now a straightforward task 
to verify that all commutation relations of Table 2 are fulfilled by the boson images 
obtained above. 

Before turning to the study of the two limits of the model, we give for completeness 
some additional information on sUg(3): 



i) The irreps of sug(3) have been studied in 
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169 


170 


171 


172 


173 


174 



T7^ , p75| , p^ , p77| 



ii) Clebsch-Gordan coefficients for sug(3) have been given in 
The Casimir operators of sUg(3) and 



L67,17e,17ah7l 



in li^. Using the Elhott quantum numbers 

A = /i - /2, 



genvalues have been given exphcitly 

183 



/2 



(27.12) 



where fi represents the number of boxes in the i-th hne of the corresponding Young 
diagram, the irreps of sUg(3) are labelled as (A,/i). The eigenvalues of the second 
order Casimir operator then read 



Co 



+ 



2A u 

h - 

3 3 



1 2 



2/i 



1 2 



while the eigenvalues of the third order Casimir operator are 
C3 = 2 



2A fi ^ 

h - + 1 

3 3 



A 2^ ^ 
- + — + 1 
3 3 



In the limit q 



1 these reproduce the ordinary results for su(3): 
2 



C2 = -(A' + /i' + A/x + 3A + 3/i), 



C, 
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(A-/i)(2A + /i + 3)(A + 2/i + 3). 



iv) A different deformation of si (3) has been studied in 



(27.13) 

(27.14) 

(27.15) 
(27.16) 



27.2. The sug(2) Urmt 

We shall study the sUq(2) limit of the model ffist, since it is technically less 
demanding. 

So far we have managed to write a boson realization of Uq(3) in terms of 3 q- 
bosons, namely ai, 02, 03. Omitting the generators involving one of the bosons, one 
is left with an sug(2) subalgebra. Omitting the generators involving 03, for example, 
one is left with A12, A21, Ni, N2, which satisfy usual sUg(2) commutation relations if 
the identifications 

J+ = /li2, J-=A2i, Jo = ^iN,-N2) (27.17) 

are made. Jq alone forms an sOg(2) subalgebra. Therefore the relevant chain of 
subalgebras is 

su,(3) D su,(2) D sOg(2). (27.18) 

The second order Casimir operator of sug(2) has been given in eq. (14.9). Substi- 
tuting the above expressions for the generators one finds 

C^2(su,(2)) = 



^Vi + I\2 



+ 1 



(27.19) 



All of the above equations go to their classical counterparts by allowing q—^ 1 



for which [x\ — > x, i.e. g-numbers become usual numbers. In the classical case out 
of the three bosons (ao, a+, a J) forming su(3), one chooses to leave out ao, the boson 
with zero angular momentum, in order to be left with the su(2) subalgebra formed 
by a_|_ and a_, the two bosons of angular momentum two. The choice of the sUq(2) 
subalgebra made above is then consistent with the following correspondence between 
classical bosons and g-bosons 

a_|_ ai, a_ a2, ao a^. (27.20) 

(We have opted in using different indices for usual bosons and g-bosons in order to 
avoid confusion.) 

In the classical case the states of the system are characterized by the quantum 
numbers characterizing the irreducible representations (irreps) of the algebras ap- 
pearing in the classical counterpart of the chain of eq. (27.18). For su(3) the total 
number of bosons N is used. For su(2) and so(2) one can use the eigenvalues of 
and Jo, or, equivalently, the eigenvalues of a^f.a_|_ + ala_ and L3 = 4Jo, for which we 
use the symbols Ud (the number of bosons with ^ng ular momentum 2) and M. Then 
the basis in the classical case can be written as 



\N,nd,M>= --^^±^^ rr^i^ r^|0 > • (27.21) 

' ' {N -ndy.{nd/2 + M/4)\{nd/2-M/A)V ^ ^ 

In the quantum case for each oscillator one defines the basis as in sec. 10. Then the 
full basis in the g-deformed case is 

(at)^-"'* (al)'^^/2+M/4 ^J^n,/2-M/4 

where N = Ni + N2 + is the total number of bosons, nd = Ni + N2 is the number 
of bosons with angular momentum 2, and M is the eigenvalue of L = AJq. Ud takes 
values from up to N, while for a given value of Ud, M takes the values ±2^^, 
±2(nrf — 2), . . . , ±2 or 0, depending on whether Ud is odd or even. In this basis the 
eigenvalues of the second order Casimir operator of sug(2) are then 



C2(su,(2))|iV, nd,M>, 



Ud 

2 



— + 1 
2 



\N,nd,M>g. (27.23) 



In the case of = 5 one can easily see that the spectrum will be composed by the 
ground state band, consisting of states with M = 0, 2, 4, 6, 8, 10 and Ud = M/2, the 
first excited band with states characterized by M = 0, 2, 4, 6 and Ud = M/2 + 2, and 
the second excited band, containing states with M = 0, 2 and nd = M/2 + 4. 

In the case that the Hamiltonian has the suq(2) dynamical symmetry, it can be 
written in terms of the Casimir operators of the chain (27.18). Then one has 



H = Eo + AC2(su,(2)) + 5C2(so,(2)), 



(27.24) 



where Eq, A, B are constants. Its eigenvalues are 



E = Eq + A 



rid 



nd 

— + 1 
2 



+ BM^. (27.25) 
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Realistic nuclear spectra are characterized by strong electric quadrupole transi- 
tions among the levels of the same band, as well as by interband transitions. In the 
mework of the present toy model one can define, by analogy to the classical case 
, quadrupole transition operators 

(5+ = difls + ala2, Q- = 02^3 + alai. (27.26) 

In order to calculate transition matrix elements of these operators one only needs 
eqs. (10.7), (10.8), i.e. the action of the g-boson operators on the g-deformed basis. 
The selection rules, as in the classical case, are AM = ±2, And = il? while the 
corresponding matrix elements are 



<N,nd+l,M±2\Q±\N,nd,M>g= W [iV - n,] [y ± ^ + 1], (27.27) 



,<N,nd-l,M± 2\Q±\N, rid, M >,= ^ [N - rid + Ijl^ T (27.28) 

From these equations it is clear that both intraband and interband transitions are 
possible. 

In order to get a feeling of the qualitative changes in the spectrum and the tran- 
sition matrix elements resulting from the g-deformation of the model, a simple cal- 
culation for a system of 20 bosons (A^ = 20) has been performed in E^. Two cases 
are distinguished: i) q real (g = e"^, with r real), ii) g a phase factor (g = e*"^, with r 
real). The main conclusions are: 

i) When g is real the spectrum is increasing more rapidly than in the classical 
case, while when g is a phase the spectrum increases more slowly than in the classical 
case, in agreement with the findings of the g-rotator model. 

ii) The transition matrix elements in the case that g is real increase more rapidly 
than in the classical case, while they increase less rapidly than in the classical case 
when g is a phase. 

iii) Transition matrix elements are much more sensitive to g-deformation than 
energy spectra. This is an interesting feature, showing that g-deformed algebraic 
models can be much more flexible in the description of transition probabilities than 
their classical counterparts. 



27.3. The S0g(3) Umit 

The classical su(3) toy model has, in addition to the above mentioned su(2) 
chain of subalgebras, an so(3) chain. However, the problem of constructing the 
sUg(3)DsOg(3) decomposition is a very difficult one. Since this decomposition is needed 



in mriRtnirt iriFr th e a -H p formerl versions of several collective models, including the El- 
liott model pi y i84P3 y i86| g^j^^ su(3) limit of the IBM, we report here the state of 
the art in this problem: 

i) As far as the completely symrne tric irreps of sUg(3) are concerned, the problem 
has been solved by Van der Jeugt E^' PT ' P9 y i90 |j i9i [ 'pj-^jg gufgces for our needs in the 
framework of the toy IBM, since only completely symmetric sUg(3) irreps occur in it. 

ii) Sciarrinoli23 started from sOg(3) and obtained a deformed gl(3) containing sOg(3) 
as a subalgebra. However, it was not clear how to impose the Hopf structure on this 
larger algebra. Trying the other way around, he found that by starting from a gig (3) 
algebra, which already possesses the Hopf structure, one loses the Hopf structure of 
the principal 3-dim subalgebra, which should have been sOq(3). 



iii) Pan 11221 and Del Sol Mesa et al t2j attacked the problem through the use of 



g-deforming functionals of sees 10, 14. 

iv) Quesne EhI started with g-bosonic operators transforming as vectors under 
sOg(3) and constructed a g-deformed u(3) by tensor coupling. 

v) g-deformed subalgebras of several g-deformed algebras have recently been stud- 
ied by Sciarrino 

vi ) A simplified version of the sOq(3) subalgebra of Uq(3) has been constructed 
in 1123. Furthermore, explicit expressions for the irreducible vector and quadrupole 
tensor operators under sOg(3) have been given and the matrix elements of the latter 
have been calculated. 



In what follows, it suffices to use the solution given in t^, since in the model under 
study only completely symmetric irreps of Uq(3) enter. Using the notation 

a+ ai, a_ 02, cts? (27.29) 

the basis states are of the form 

n+non_ >= ^ + ^ ^> >, 27.30 

W]\[n,]\[n_]\ 



with ai|0 >= and A^j|n+non_ >= nj|n+non_ >, where i = +,0,—. The principal 
subalgebra sOg(3) is then generated by 

Lo = N+-N_, (27.31) 



L+ = g^— l^«yg^+ + g-^+ a^ao + g^+-|^-' Jg^- + g-^- , (27.32) 



L_ = 4a+g^— ^^«V9^v^Tg^ + g^+-^^°VV^^+9^«-«o, (27.33) 
satisfying the commutation relations 

[Lo, L±] = ±L±, L_] = [2Lo]. (27.34) 

Lq alone generates then the sOg(2) subalgebra. Therefore the relevant chain of subal- 
gebras is 

sug(3) D sOg(3) D so5(2). (27.35) 



The sOg(3) basis vectors can be written in terms of the vectors of eq. (27.30) as 



\v{N,L,M) >= ^-[{^+M){^+Af-l)]/4^ 



[A^ + L]!![2L + 1][L + M]\[L - M]\ 
\ [N - L\\\[N + L + l\\ 



(2L-i)x/2 g{N-L)/2 \x,L + M - 2x,x - M > (27 36) 

'[2x]\\[L + M - 2x]\[2x - 2M]U' 



where 



s = (aS)^g^-+^-+^ 



[2iV+][2iV_j t t^-^o-i 



\ [NAN-] 



a\alq (27.37) 



X takes values from max(0,M) to [{L + M)/2] in steps of 1, L = iV, iV - 2, . . . , 1 
or 0, M = -L, -L + 1, . . . , +L, and [2x\\\ = [2x][2x -2]... [2]. The action of the 
generators of sOq(3) on these states is given by 

Lo\v{N, L, M) >= M \v{N, L, M) >, (27.38) 



L± \v{N,L,M) >= ^J[Lt M][L±M + 1] \v{N, L, M ±1) > . (27.39) 
The second order Casimir operator of sOg(3) has the form 

C2(so,(3)) =L^ = L_L+ + [Lo][Lo + 1]. (27.40) 

Its eigenvalues in the above basis are given by 

(:72(so,(3)) \v{N, L, M) >= [L] [L + 1] \v{N, L,M)> . {27 Al) 

All of the above equations go to their classical counterparts by allowing q—^ 1, 
for which [x] — > x, i.e. g-numbers become usual numbers. In the classical case the 
states of the system are characterized by the quantum numbers characterizing the 
irreducible representations (irreps) of the algebras appearing in the classical counter- 
part of the chain of eq. (27.35). For su(3) the total nuniber of bosons N is used. For 
so(3) and so(2) one can use L and M, respectively. In ll^, however, the eigenvalue of 
L'q = 2Lq is used, which is M' = 2M. 

Since the rules for the decomposition of the totally symmetric Uq(3) irreps into 
sOg(3) irreps are the same as in the classical case, it is easy to verify that for a system 
with = 6 the spectrum will be composed by the ground state band, consisting 
of states with M' = 0, 2, 4, 6, 8, 10, 12 and L = N, the first excited band with 
states characterized by M' = 0, 2, 4, 6, 8 and L = N — 2, the second excited band, 
containing states with M = 0, 2, 4 and L = A^ — 4, and the third excited band, 
containing a state with M' = and L = N — 6. 

In the case that the Hamiltonian has the sOq(3) dynamical symmetry, it can be 
written in terms of the Casimir operators of the chain (27.35). Then one has 



H = Eo + AC2{sOg{3)) + 5C2(so,(2)), 



(27.42) 



where Eq, A, B are constants. Its eigenvalues are 



E = Eo + A[L][L + 1] + BM^ 



(27.43) 



It is then clear that in this simple model the internal structure of the rotational bands 
is not influenced by g-deformation. What is changed is the position of the bandheads. 

We turn now to the study of electromagnetic transitions. In the present limit 
one can define, by analogy to the classical case, quadrupole transition operators Q± 
proportional to the sOg(3) generators L± 



(27.44) 



In order to calculate transition matrix elements of these operators one only needs eq. 
(27.39). The selection rules, as in the classical case, are AM' = ±2, AL = 0, i.e. 
only intraband transitions are allowed. The relevant matrix elements are 



< v{N, L, M' + 2)\Q+\v{N, L, M') >-- 



~ ~2~. 



■ M' ■ 



(27.45) 



< v{N, L, M' - 2)\Q_\v{N, L, M') >-- 



■ M' 



(27.46) 



In order to get a feeling of the qualitative changes in the spectrum and the tran- 
sition matrix elements resulting fro rn th e g-deformation of this limit of the model, a 
simple calculation has been done in 1123. Again the cases of q being real or q being a 
phase factor have been considered. The main conclusions are: 

i) g-deformation influences only the position of bandheads, while it leaves the 
internal structure of the bands intact. 

ii) When q is real the bandheads are increasing more rapidly than in the classical 
case, while when g is a phase the bandheads increase more slowly than in the classical 
case. This result is in qualitative agreement with the findings of the suq(2) model. 

iii) Transition matrix elements in the case that g is real have values higher than 
in the classical case, while they have values lower than in the classical case when g is 
a phase. 



28. The question of complete breaking of symmetries and some applica- 
tions 

In the cases examined so far the g-deformed symmetries considered were close 
to their classical counterparts, to which they reduce for r = (g = 1), since the 
values of r were relatively small. One can then argue that results similar to the 
ones provided by the quantum symmetries can also be obtained from the usual Lie 
symmetries through the addition of suitable perturbations. What can be very useful 
is to start with one limiting symmetry and, through large deformations, reach another 
limiting symmetry. We shall refer to this as the complete breaking of the symmetry. 



In this way one could hope to "bridge" different Lie symmetries through the use 
of g-deformations, providing in addition new symmetries in the regions intermediate 
between the existing Lie ones. 

The question of complete breaking of the^ymmetry in the framework of the toy 



under study has been studied by Cseh Gupta and Del Sol Mesa et al. 



Cseh t£j started with the sug(2) (vibrational) limit (in a form different from the 
one used above) and tried to reach the sOq(3) (rotational) limit. He noted that for q 
being a phase factor this is not possible, while for real q some rotational features are 
obtained, but without all of the requirements for rotational behaviour being satisfied 
simult aneously. 

Gupta started with the suq(2) limit, in the form given above. He noted that 
for q being a phase factor a recovery of the su(3) symmetry occurs (see also the next 
paragraph), while for real g, and even better for q complex (g = with s = a + ih), 
the sOg(3) hmit is indeed reached. 



Del Sol Mesa et al. E£3 considered the 0^(8) limit of the model, since it corresponds 
to the symmetry of a g-deformed version of the spherical Nilsson Hamiltonian with 
spin-orbit coupling term. They found that for g being a phase factor (g = e*^) and 
for r obtaining values in the region 0.5 < r < 2 the u(3) symmetry, which is broken 
in the initial model because of the presence of the spin-orbit term, is recovered. This 
offers a way of recoverin g the u(8) symmetry alternative to the o ne d eveloped for the 



spherical Nilsson model r^^ H ^"^ and the deformed Nilsson model EiS'E!^ through the use 



of appropriate unitary operators. 

Complex deformations have also been used in c2a in the framework of a deformed 
u(2) model, possessing the u(2)du(1) and u(2)do(2) chains. Again it has been found 
that complex deformations can bridge the two limiting symmetries. 

Possible complete breaking of the symmetry has also been studied in the frame- 
work of the g-deformed version of the full Interacting Boson Model (see sec. 29). 

A problem associated with complex deformations as the ones considered above is 
that the energy eigenvalues become complex as well. A way to avoid this problem 
has been introduced recently by Jannussis and collaborators E^IE^ 
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Finally, the Oq(3) limit of the model has been used for describing the ^^O + a 
cluster states in ^°Ne It turns out that an improved description of the energy 
spectrum and the a-particle spectroscopic factors occurs for g = e°'^^^. 

29. g-deformation of the Interacting Boson Model (IBM) 



The Interacting Boson Model E^f 5ti^|i57| ^g^g |i2y,|i3S||'Qj, rg^gnt overviews) is the most 
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138 



popular algebraic model of nuclear structure. It describes the collective properties 
of medium-mass and heavy nuclei away from closed shells in terms of bosons, which 
correspond to correlated valence fermion pairs. In its simplest form, called IBM-1, 
only s ( J = 0) and d {J = 2) bosons are used. The overall symmetry of the model is 
u(6), possessing three limiting symmetries: the u(5) (vibrational) limit, corresponding 



to the chain of subalgebras 



u(6) D u(5) D o(5) D o(3), 
u(6) D su(3) D o(3), 



u(6) D 0(6) D 0(5) D 0(3). 



(29.1) 
(29.2) 
(29.3) 



the su(3) (rotational) hmit, characterized by the chain 



and the 0(6) (7-unstable) hmit, for which the relevant chain is 



If one of these dynamical symmetries is present, the Hamiltonian can be written in 
terms of the Casimir operators of the algebras appearing in the relevant chain. Thus 
the Hamiltonian can be analytically diagonalized in the corresponding basis. This 
is a great advantage of IBM and its numerous generalizations: they provide us with 
a large number of exactly soluble models, the predictions of which can be directly 
compared to experiment, without any need for lengthy numerical calculations. 

From what we have already seen in sec. 27, it is worth examining if a g-deformed 
version of the IBM has any advantages in comparison to the standard version. In 
order to accomplish this, one has to construct the g-analogues of the three chains 
mentioned above. The difficulties associated with the su(3) chain have already been 
discussed in subsec. 27.3. In what follows we are going to focus attention on the 0(6) 
chain, for which the relevant construction has been carried out 0. The technique 
used is based on the notion of complementary subalgebras, which is explained in 
detail in Ei^, while here only final results will be reported. We only mention here that 
the notion of complementary subalgebras was introduced by Moshinsky and Quesne 
^lU'ElB'tij. Two subalgebras Ai and A2 of a larger algebra A are complementary within 
a definite irrep of A, if there is an one-to-one correspondence between all the irreps 
of Ai and A2 contained in this irrep of A. 

In the 0(6) limit of IBM the Hamiltonian is 

H = Eo + /?C2(o(5)) + 7^^2(0(3)) + r/C2(o(6)). (29.4) 

The eigenvalues of the energy in the relevant basis have already been given in subsec. 
26.1. 

Using the notion of complementarity it turns out that, instead of using the 0(6) 
chain mentioned above, it suffices to study the chain 

su"'^(l, 1) (g) so(6) D su'^(l, 1) ® so(5) D so(3), (29.5) 

where su*'^(l,l) is the algebra closed by the pair operators formed out of the s and 
d bosons, while su'^(l,l) is the algebra closed by the pair operators formed out of 
d bosons alone. (Details on the basis for symmetric irreps of su(l,l)®o(5) can be 
found in 0.) The irreps of su*'^(l,l) are characterized by the same quantum numbers 
as the irreps of 0(6) in the 0(6) chain of the IBM, while the irreps of su'^(l,l) are 



characterized by the same quantum numbers as the irreps of o(5) in the 0(6) hmit of 
IBM. Therefore in the Hamihonian one can use the Casimir operators of the su'"^(l,l), 
su'^(l,l) and su(2) subalgebras (the deformed versions of which are well known, as 
seen in sees 14-16) instead of the Casimir operators of 0(6), o(5), o(3) respectively. 
Keeping the same notation as in eq. (26.3) the final result reads 



E(Ar, d, r, i/A, J, Mj) = Eo + /58 



+72[J]JJ + l], + r/8 
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(29.6) 



where the free parameters have been chosen so that the present equation reduces to 
its classical counterpart for g — ^ 1. For the ground state band then the analog of eq. 
(26.4) is 



E{J) = Eo + /3'8[J],i/4[J + 6],i/4 + 72[J] JJ + l]c + v'8[N]M^ + 4]< 



1/2, 



where the identities 



X 



-1/2N 



(29.7) 



(29. 



(29.9) 



have been used and jS', rj' are related to (3, rj and q in an obvious way. We remark that 
the Casimir operator of sUg(l,l), which is complementary to o(5) in the undeformed 
case, leads to a term of the form [Jjq'f'^ + with q' = q 
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the question has been studied if large values of the deformation 
parameter can lead us from the 0(6) limit to the su(3) (rotational) or u(5) (vibrational) 
limits, so that complete breaking of the symmetry, in the sense of sec. 28, could be 
obtained. It turns out that for q real the spectrum of the ground state band goes 
towards the rotational limit, while q being a phase factor leads towards the vibrational 
limit. Many more detailed studies, of both spectra and electromagnetic transition 
probabilities, are required before such a claim can be made. 

A different method for constructing th e g- deformed versions of the u(5) and 0(6) 
limits of the IBM has been used by Pan £13. The method is based on the use of 
g-deforming functionals (see sees 10, 14). The same method has been used in 1123 for 
studying the g-deformed version of the su(3)Dso(3) decomposition. The final result 
for the energy eigenvalues in the 0(6) case is similar to the one reported in eq. (29.6), 
the main difference being that different deformation parameters are allowed in each 
of the three deformed terms in the rhs. Some comparisons of the model predictions 
for spectra and B(E2) values to the experimental data have been performed il3. 

It is clear that several deformed versions of the IBM can be constructed, providing 
us with a large number of exactly soluble models. In order to demonstrate their 
usefulness, one has to show that by deforming the model one gets some advantages 
over the classical (non-deformed) version. One way to achieve this is the use of 



paxamfiter-independent tests based on systematics of the data, hke the ones used in 
El3'E13'E23 for the usual IBM. It is also desirable for the deformation parameter to be 
associated with some physical quantity, as in the case of the sUg(2) model. Much 
work is still required in these directions. Some mathematical results which can be 
useful in these efforts are reported below: 

i) Casimir operators f or !■ 
of sOq(5) can be found in 



Uq(n) have been given in while the quadratic Casimir 
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ii) Raising and lowering operators for Ug(n) have been given by Quesne 

iii) Irreps of Ug(m+n) in the Uq(m)©Ug(n) basis have been constructed inE2J, while 
generalized g-bosonic operators acting in a tensor product of m Fock space s have been 
constructed as double irreducible tensors with respect to Ug(m)©Ug(n) in "^^^ 



30. Deformed versions of other collective models 



The Moszkowski model 



is a schematic two-level model which provides a de- 
scription of the phase transition from the vibrational regirpe to the rotational one. 



A g-deformed version of the model has been developed cHa and the RPA modes in 
it have been discussed iH. Furthermore, the g-deformed Moszkowski model with 
cranking has been studied in the mean field approxirnation and the relation between 
g-deformation and temperature has been discussed It should be noticed here 
that quantum algebraic techniques have also been found useful in describing thermal 
effects in the framework of the g-deformed Thouless model for superconductivity c^^ 



The Lipkin-Meshkov-Glick (LMG) model is an exactly soluble schematic shell 
model, g-deformed versions of the 2- level LMG model (in terms of an sUg(2) algebra) 
and of the 3- level LMG model (in terms of an sUg(3) algebra) E3 have 
been developed. 



31. Fermion pairs as deformed bosons: approximate mapping 

We have seen so far that several quantum algebraic phenomenological models 
have been proposed for the description of nuclear collective properties. These models 
make use of g-deformed bosons, which satisfy commutation relations differing from 
the standard boson commutation relations, to which they reduce in the limit g — 1. 

On the other hand, it is known that vibrational nuclear spectra, which are de- 
scribed in the simplest way by a pairing Hamiltonian, show anharmonicities (see also 
sec. 26), described, for example, by the Anharmonic Vibrator Model (AVM) llH 
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(31.1) 



which can be extended 



E{J) = aJ + hJ{J -2). 
In the framework of the single-j shell 



iel 



to several non-degenerate j-shells Eiil'cij'Bfa these anharmonicities are related to the 
fact that correlated fermion pairs satisfy commutation relations which resemble boson 
commutation relations but in addition include corrections due to the presence of the 
Pauli principle. This fact has been the cause for the development of boson mapping 



techniques (see the recent reviews by Klein and Marshalek and Hecht ll^ and ref- 
erences therein), by which the description of systems of fermions in terms of bosons 
is achieved. In recent years boson mappings have attracted additional attention in 
nuclear physics as a necessary tool in providing a theoretical justification for the suc- 
cess of the phenomenological Interacting Boson Model and its various extentions, 
in which low lying collective states of medium and heavy mass nuclei are described 
in terms of bosons. 

From the above observations it is clear that both q-bosons and correlated fermion 
pairs satisfy commutation relations which resemble the standard boson commutation 
relations but they deviate from them, due to the g-deformation in the former case 
and to the Pauli principle in the latter. A question is thus created: Are g-bosons 
suitable for the approximate description of correlated fermion pairs? In particular, 
is it possible to construct a boson mapping in which correlated fermion pairs are 
mapped onto g-bosons, in a way that the g-boson operators approximately satisfy the 
same commutation relations as the correlated fermion pair operators? In this section 
we show for the simple case of su(2) that such a mapping is indeed possible. 

31.1. The single- j shell model 
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pair and multipole operators as 

4m = 7f E {jmjm'\JM)a]^a]^,, (31.2) 

V ^ mm' 

BJM = ^=i=T E U^J - mVM)(-l)^-'a]„a,w, (31.3) 

V ^'J + ^ mm' 

with the following definitions 

Ajm = [AIm]\ bI,, = [Bjm]^. (31.4) 

In the above {a,jm) are fermion creation (annihilation) operators and {jmjm'lJM) 
are the usual Clebsch-Gordan coefficients. 

The pair and multipole operators given above satisfy the following commutation 
relations: 

[4m,4'M'] = 0, (31.5) 



[AjM, 4'M'] = hj'hiM' - 2 ^(-1)2^+^V(2J + 1)(2 J' + 1)(2J" + 1) 

J" 

{J - MJ'M'\J"M' - M) S^"^^ ^' ^"|5j»,M'-M, (31.6) 

[BjM, 4'M'] = E 2V2J' + - MJ'M'\J"M' - M) 

J" 



j j j j^J",M'-M 2 ' W^-') 

[BjM,Bj,M'] = E(-i)''"'"[i - {-iy^'''-'"]V2J^ 

J" 

( JMJ'MV"M + M') I ^ ^ '^^^Bjn^M+M', (31.8) 

in which the curly brackets are the usual 6-j symbols. These are the commutation 
relations of the so(2(2j+l)) algebra. 

31.2. Fermion pairs of zero angular momentum 

In the present subsection we will restrict ourselves to fermion pairs coupled to 
angular momentum zero. The relevant commutation relations take the form 

[^0,4] = 1-^, [^,4] = 4, [^,Ao] = -Ao, (31.9) 
where Np is the number of fermions, 2Q = 2j + 1 is the size of the shell, and 

Bo = NF/V2n. (31.10) 

With the identifications 

j+ = VnAl j_ = VnAo, Jo = (31.11) 

eqs (31.9) take the form of the usual su(2) commutation relations 

[J+,J_]=2Jo, [Jo,J+] = J+, [Jo,J-] = -J- (31.12) 
An exact boson mapping of the su(2) algebra is given in 



4 = 4yl-^, Ao = yi-^ao, Np = 2no, (31.13) 

where aj (ao) are boson creation (annihilation) operators carrying angular momentum 
zero and no is the number of these bosons. 

The simplest pairing Hamiltonian one can consider has the form 

H = -GnAlAo. (31.14) 

The Casimir operator of su(2) can be written as 



while the pairing energy takes the form 



E _Nf Nl ^ Nf 
(-GO) "2 4Q + 2Q- ^ 

Our aim is to check if there is a boson mapping for the operators Aq and Np 
in terms of g-deformed bosons, having the foUowing properties: 

i) The mapping is simpler than the one of eq. (31.13), i.e. to each fermion 
pair operator ^4^, Aq corresponds a bare g-boson operator and not a boson operator 
accompanied by a square root (the Pauli reduction factor). 

ii) The commutation relations (31.9) are satisfied up to a certain order. 

ii) The pairing energies of eq. (31.16) are reproduced up to the same order. 

31.3. The q-deformed oscillator: A story of failure 

In the case of the g-deformed harmonic oscillator (sec. 10), the commutation 
relation 

[a,a^] = [N + 1]-[N] (31.17) 
for q being a phase can be written as 

(2jV+l)r 

[a,a^]- I ■ (31.18) 

cos I 

In physical situations r is expected to be small (i.e. of the order of 0.01). Therefore 
in eq. (31.18) one can take Taylor expansions of the functions appearing there and 
thus find an expansion of the form 

2 4 

[a, a^] = 1 - ^{N"^ + N) + ^{N^ + 2N'^ - N) - . . . . (31.19) 

We remark that the first order corrections contain not only a term proportional to N , 
but in addition a term proportional to A?"^, which is larger than N . Thus one cannot 
make the simple mapping 

Ao^a, Al^a^, Np ^ 2N, (31.20) 

because then one cannot get the first of the commutation relations (31.9) correctly 
up to the first order of the corrections. The same problem appears in the case that q 
is real as well. In addition, by making the simple mapping of eq. (31.20) the pairing 
Hamiltonian can be written as 



(31.21) 



In the case of small r, one can take Taylor expansions of the functions appearing 
in the definition of the gr-numbers (eq. (2.2) or (2.3)) and thus obtain the following 
expansion 

2 4 

\N]^ N±—(N - N^) + —(7N - lON^ + 3N^) 
^ ^ 6 ^ ^ 360^ ^ 

± (31iV - 49iV^ + 21N^ - 3iV^) + . . . , (31.22) 

15120 

where the upper (lower) sign corresponds to q being a phase factor (real). We remark 
that while the first order corrections in eq. (31.16) are proportional to Np and Np, 
here the first order corrections are proportional to N and A^^. Thus neither the 
pairing energies can be reproduced correctly by this mapping. 

31.4- The Q-deformed oscillator: A story of success 

In the case of the Q-oscillator of sec. 11, however, the commutation relation 
among the bosons is 

[6,6t] = g^. (31.23) 
Defining Q — this can be written as 

= 1 + TiV+ ^— + — — + . . . . (31.24) 
2 

We remark that the first order correction is proportional to N. Thus, by making the 
boson mapping 

Al^h\ Ao^b, Nf^2N, (31.25) 

one can satisfy the first commutation relation of eq. (31.9) up to the first order of 
the corrections by determining T — —2/Q. 

We should now check if the pairing energies (eq. (31.16)) can be found correctly 
up to the same order of approximation when this mapping is employed. The pairing 
Hamiltonian in this case takes the form 

= fttj, ^ [^]^. (31.26) 

Defining Q = it is instructive to construct the expansion of the Q-number of 
eq. (6.1) in powers of T. Assuming that T is small and taking Taylor expansions in 
eq. (6.1) one finally has 

[N]q = N+-{N'^ -N) + —{2N^ - 3N^ + 1) + —{N^- 2N^ + N^) + . . . (31.27) 

Using the value of the deformation parameter T — —2/fl, determined above from the 
requirement that the commutation relations are satisfied up to first order corrections, 
the pairing energies become 

E N^-N 2Ar3_3Ar2 + i n'^-2N^ + N^ 

+ — + .... (31.28) 



-on n 3f72 3j^3 



The first two terms in the rhs of eq. (31.28), which correspond to the leading term plus 
the first order corrections, are exactly equal to the pairing energies of eq. (31.16), 
since Np — >■ 2N. We therefore conclude that through the boson mapping of eq. 
(31.25) one can both satisfy the fermion pair commutation relations of eq. (31.9) and 
reproduce the pairing energies of eq. (31.16) up to the first order corrections. 
The following comments are also in place: 

i) By studying the spectra of the two versions of the g-deformed harmonic oscilla- 
tor, given in eqs. (10.10) and (11.10), one can easily draw the following conclusions: 
when compared to the usual oscillator spectrum, which is equidistant, the spectrum 
of the g-oscillator is getting shrunk for q being a phase, while the spectrum of the Q- 
oscillator gets shrunk when T < 0. In a similar way, the spectrum of the g-oscillator 
gets expanded for q real, while the spectrum of the Q-oscillator gets expanded for 
T > 0. In physical situations (sees 19-23) it has been found that the physically in- 
teresting results are gotten with q being a phase. Thus in the case of the Q-oscillator 
it is the case T < the one which corresponds to the physically interesting case. As 
we have already seen, it is exactly for T = —2/Q < that the present mapping gives 
the fermion pair results. 

ii) It should be recalled that the pairing model under discussion is studied under 
the assumptions that the degeneracy of the shell is large {Q » 1), that the number 
of particles is large (A^ >> 1), and that one stays away from the center of the shell 
[Q — N = 0{N)). The accuracy of the present mapping in reproducing the pairing 
energies has been checked in cia, where results for f2 = 11 (the size of the nuclear fp 
major shell), Q = 16 (the size of the nuclear sdg major shell) and Q = 22 (the size of 
the nuclear pfh major shell) are reported, along with results for the case = 50 (as 
an example of a large shell). In all cases good agreement between the classical pairing 
model results and the Q-Hamiltonian of eq. (31.26) is obtained up to the point at 
which about 1/4 of the shell is filled. The deviations observed near the middle of the 
shell are expected, since there the expansion used breaks down. 

We have thus shown that an approximate mapping of the fermion pairs coupled 
to angular momentum zero in a single-j shell onto suitably defined g-bosons (the Q- 
bosons) is possible. The su(2) commutation relations are satisfied up to the first order 
corrections, while at the same time the eigenvalues of a simple pairing Hamiltonian 
are correctly reproduced up to the same order. The small parameter of the expansion, 
which is T (where Q = e^), turns out to be negative and inversely proportional to 
the size of the shell. 

The present results are an indication that suitably defined g-bosons could be used 
for approximately describing systems of correlated fermions under certain conditions 
in a simplified way. The construction of g-bosons which would exactly satisfy the 
fermion pair su(2) commutation relations will be undertaken in the following section. 

32. Fermion pairs as deformed bosons: exact mapping 

From the contents of the previous section, the following question is created: Is it 



possible to construct a generalized deformed oscillator (as in sec. 12) using deformed 
bosons in such a way that the spectrum of the oscillator will exactly correspond 
to the pairing energy in the single-j shell model, while the commutation relations 
of the deformed bosons will exactly correspond to the commutation relations of the 
correlated fermion pairs in the single-j shell under discussion? In this section we show 
that such an oscillator can indeed be constructed cl3 by using the method of sec. 12. 

32.1. An appropriate generalized deformed oscillator 

We apply the procedure of sec. 12 in the case of the pairing in a single-j shell 
mentioned before. The boson number is half the fermion number, i.e. = Np/2. 
Then eq. (31.16) can be written as 

E N'^ N , , 

^-T^+o- (32.1) 



One can then use a generalized deformed oscillator with structure function 

F{N) = a^a = N- — + -. (32.2) 

In addition one has 

P^r,^r,^aa^^N+l-^-^ + ^. (32,3) 

What we have constructed is a boson mapping for the operators Aq, Aq, Np: 

Ao^a, Al^a\ Np ^ 2N. (32.4) 

From eq. (32.2) it is clear that this mapping gives the correct pairing energy. In 
addition one has 

2N 

[a, at] = F{N + 1) - F{N) = 1 - (32.5) 

in agreement to the first commutation relation of eq. (31.9). Thus the correct com- 
mutation relations are also obeyed. (The last two commutation relations of eq. (31.9) 
are satisfied because of (12.1).) 

As we have already seen in the previous section, an exact hermitian boson mapping 
for the su(2) algebra is known to have the form of eq. (31.13). In this mapping the 
Pauli principle effects are carried by the square roots accompanying the ordinary 
boson operators, while in the mapping of eq. (32.4) the Pauli principle effects are 
"built in" the deformed bosons. 

The generalized oscillator obtained here has energy spectrum 

En = liFiN) + F(iV + 1)) = iV + i - ^, (32.6) 



which is the spectrum of an anharmonic oscillator. 



32.2. Related potentials 

The classical potential giving the same spectrum, up to first order perturbation 
theory, can be easily determined (see also subsecs 13.1, 13.2). The potential 

V{x) = Kx'^ + Xx\ (32.7) 

is known to give in first order perturbation theory the spectrum 

En = K{2n + 1) + A(6n^ + 6n + 3) = {2k + 6A)(n +h + GAn^. (32.8) 
Comparing eqs. (32.6) and (32.8) one finds 

Then the classical potential giving the same spectrum, up to first order perturbation 
theory, as the generalized oscillator determined here, is 

It is therefore demonstrated that the Pauli principle effects in a single-j shell with 
pairing interaction are equivalent to an anharmonicity. 

The generalization of the results obtained in this section for the pairing hamilto- 
nian to any anharmonic oscillator is straightforward. For example, the potential 

y(a;) = + Ax^ + /xx^ + (32.11) 

is known to give up to first order perturbation theory the spectrum 

En = K{2n + 1) + A(6n2 + 6n + 3) 

+//(20n^ + 30n^ + 40n + 15) + C(70n^ + 140n^ + 350/1^ + 280n + 105), (32.12) 
which can be rewritten in the form 

En={n+{n+ 1)){k + 5/x) + (n' + (n + 1)')(3A + 70^ 

+(71^ + (n + 1)^)(10a() + {n^ + (n + 1)^)(350. (32.13) 
Taking into account eq. (12.11), from eq. (32.13) one gets 

= (k + 5//)n + (3A + 700n^ + (10/i)n^ + (35^^^ (32.14) 

For /i = ^ = and A given from eq. (32.9), the results for the pairing problem are 
regained. 



It is worth mentioning at this point that the energy spectrum of the generahzed 
oscillator corresponding to the pairing correlations (eq. (32.6)) can be rewritten as 




On the other hand, it is known that for the modified Poschl-Teller potential (see also 
subsec. 13.2) 

V{x) = Dtanh'^{x/R), (32.16) 

25^ 



the energy spectrum is given by 

It is thus clear that the energy spectrum of the generalized oscillator studied here 
can be obtained from the modified Poschl-Teller potential for special values of the 
potential depth D. 

It is also worth remarking that the "structure function" F{N) of the generalized 
oscillator obtained here (eq. (32.2)) can be written as 

N 

F{N) = —{n + l-N), (32.18) 

which is similar to the one of the para-fermionic oscillator of Ohnuki and Kamefuchi 
@ (see also sees 12, 18). 

In summary, we have constructed a generalized deformed oscillator which satis- 
fies the same commutation relations as fermion pair and multipole operators of zero 
angular momentum in a single-j shell, and, in addition, reproduces the pairing en- 
ergy of this shell exactly. We have thus demonstrated that an exact hermitian boson 
mapping of a system of angular-momentum-zero fermion pairs in terms of bare de- 
formed bosons can be constructed, while in the usual case the ordinary bosons are 
accompanied by square roots due to the Pauli principle effects. The oscillator corre- 
sponding to the pairing problem has a spectrum which can be reproduced up to first 
order perturbation theory by a harmonic oscillator with an anharmonicity. The 
construction of a generalized deformed oscillator corresponding to any anharmonic 
oscillator has also been achieved. 

The results obtained in this section indicate that deformed bosons might be a 
convenient tool for describing systems of fermion pairs under certain conditions. The 
generalisation of the results obtained here to fermion pairs of nonzero angular momen- 
tum, which will allow for a fuller treatment of the single-j shell in terms of deformed 
bosons, is a very interesting problem. 



33. The seniority scheme 



In the previous two sections we have seen how correlated fermion pairs of zero 
angular momentum can be mapped onto deformed bosons. It is however known 
that pairs of non-zero angular momentum play an important role in the formation of 
nuclear properties. In the present section a first step in the direction of describing 
the J 7^ pairs in terms of deformed bosons is taken. 



33.1. Uncovering a dynamical symmetry 

In the usual formulation of the theory of pairing in a single-j shell i^, fermion 
pairs of angular momentum J = are created by the pair creation operators 

S'^ = X/ {~^y^^'^]m'^\-m,^ (33.1) 



r?i>0 

where are fermion creation operators and 2VL = 2j + 1 is the degeneracy of the 
shell. In addition, pairs of nonzero angular momentum are created by the — 1 
operators 

bI=Y: (-l)^+™(jmj - m| J0)4^4_^, (33.2) 

where {jmj — m| JO) are the usual Clebsch Gordan coefficients. The fermion number 
operator is defined as 

Np = ^ ^ O^jm^im — ^ ^ i^^jm^jm ~\~ dj—m^j—rn) ■ (33.3) 
m m>0 

As we have already seen, the J = pair creation and annihilation operators satisfy 
the commutation relation 

[^,^t] = i_^, (33.4) 

while the pairing Hamiltonian is 

H = -GQS^S. (33.5) 

The seniority Vp is defined as the number of fermions not coupled to J = 0. If only 
pairs of J = are present (i.e. Vp = 0), the eigenvalues of the Hamiltonian are (as 
already seen in eq. (31.16)) 

Ei^.,V..O)^-Gn[^,^-§). ,33.6) 

For non-zero seniority the eigenvalues of the Hamiltonian are 

E{Nf, Vf) = - j(A^F - VF)i2n - Np - Vp + 2). (33.7) 



We denote the operators Np, Vp and their eigenvalues by the same symbol for sim- 
plicity. 

In subscc. 31.4 it has been proved that the behaviour of the J = pairs can be 
described, up to first order corrections, in terms of Q-bosons. In particular, making 
the mapping 

S^^b\ S^b, Nf^2N, (33.8) 
the relevant pairing Hamiltonian of eq. (33.5) becomes 

H{N, V = 0) = -GQb^b - -Gn[N]Q, (33.9) 

which coincides with eq. (33.6) up to first order corrections in the small parameter, 
which is identified as T = —2/Q. Furthermore, the Q-bosons satisfy the commutation 
relation of eq. (31.24), which coincides with eq. (33.4) up to first order corrections 
in the small parameter, which is, consistently with the above finding, identified as 
T = —2/Q. Therefore the fermion pairs of J = can be approximately described as 
Q-bosons, which correctly reproduce both the pairing energies and the commutation 
relations up to first order corrections in the small parameter. 

For the case of nonzero seniority, one observes that eq. (33.7) can be written as 

(33.0) 

i.e. it can be separated into two parts, formally identical to each other. Since the 
second part (which corresponds to the J — pairs) can be adequately described 
by the Q-bosons b, b^, and their number operator N, as we have already seen, it is 
reasonable to assume that the first part can also be described in terms of some Q- 
bosons d, S, and their number operator V (with Vp 2V), satisfying commutation 
relations similar to eqs (11.2) and (11.3): 

[V, = d\ [V, d] = -d, dd^ - QSd = 1. (33.11) 

From the physical point of view this description means that a set of Q-bosons is used 
for the J = pairs and another set for the J 7^ pairs. The latter is reasonable, since 
in the context of this theory the angular momentum value of the J 7^ pairs is not 
used explicitly. The J 7^ pairs are just counted separately from the J = pairs. A 
Hamiltonian giving the same spectrum as in eq. (33.10), up to first order corrections 
in the small parameter, can then be written as 

H{N,V)^GQ{\V]q-\N]q). (33.12) 

Using eq. (31.27) it is easy to see that this expression agrees to eq. (33.10) up to first 
order corrections in the small parameter T — —2/Q. 
Two comments concerning eq. (33.12) are in place: 

i) In the classical theory states of maximum seniority (i.e. states with N = V) 
have zero energy. This is also holding for the Hamiltonian of eq. (33.12) to all orders 
in the deformation parameter. 



ii) A landmark of the classical theory is that E{N, V)—E{N, = 0) is independent 
of N. This also holds for eq. (33.12) to all orders in the deformation parameter. 

Knowing the Schwinger realization of the sUg(2) algebra in terms of g-bosons (sec. 
15), one may wonder if the operators used here close an algebra. It is easy to see that 
the operators bid, d'^b and — ^ do not close an algebra. Considering, however, the 
operators (see " with p = 1) 



J+ = b^Q-''/^d, J. = d^Q-'''\ J, = -{N-V), 

one can easily see that they satisfy the commutation relations 00 

[Jo, J±] = ±J±, J+J- - Q-'J-J+ = [2Jo]q. 
Using the transformation 

Jo = Jo, J+ = Q(^/')(^«-^/') J+, J- = j_g(i/2){jo-i/2)^ 

one goes to the usual sug(2) commutation relations 

[Jo,J±] = ±J±, [J+, J_] = [2Jo], 



(33.13) 



(33.14) 



(33.15) 



(33.16) 



where = Q. One can thus consider eq. (33.14) as a rewriting of the algebra sUg(2), 
suitable for boson realization in terms of Q-bosons. 

It is clear that + is the first order Casimir operator of the uq(2) algebra 
formed above (since it commutes with all the generators given in eq. (33.13)), while 
A^ — \^ is the first order Casimir operator of its uqII) subalgebra, which is generated 
by Jo alone. Therefore the Hamiltonian of eq. (33.12) can be expressed in terms of 
the Casimir operators of the algebras appearing in the chain uq(2) D uq(1) as 



E{N, V) = Gn{ 



Ci(uq(2))-Ci(uq(1)) 



Ci(uq(2))+Ci(uq(1)) 



(33.17) 



Q 



i.e. the Hamiltonian has a uq(2) D uq(1) dynamical symmetry. 



33.2. Comparison to experiment 

In the construction given above we have shown that Q-bosons can be used for the 
approximate description of correlated fermion pairs in a single-j shell. The results 
obtained in the Q-formalism agree to the classical (non-deformed) results up to first 
order corrections in the small parameter. However, the Q-formalism contains in 
addition higher order terms. The question is then born if these additional terms are 
useful or not. For answering this question, the simplest comparison with experimental 
data which can be made concerns the clas sic exa mple of the neutron pair separation 
energies of the Sn isotopes, used by Talmi B3'E3. 



In Talmi's formulation of the pairing theory, the energy of the states with zero 
seniority is given by 

E{N\i = Wo + ^^^^~^^ A, (33.18) 

where N is the number of fermion pairs and Vq, A are constants. We remark that 
this expression is the same as the one in eq. (33.6), with the identifications 

A/(2Vo) = A = 2G, Nf = 2N. (33.19) 

The neutron pair separation energies are given by 

AE{N + = E{N + - E{N),i = Vo[l + yN^ . (33.20) 

Thus the neutron pair separation energies are expected to decrease hnearly with 
increasing N . (Notice from eq. (33.19) that A/Vq < 0, since f2 > 0.) A similar linear 
decrease is predicted also by the Interacting Boson Model EhI. 

In our formalism the neutron pair separation energies are given by 

AE{N + 1)q = -Gfi([iV + 1]q - [N]q) = -Gnq^ = -Gfie™. (33.21) 

Since, as we have seen, T is expected to be —2/Q, i.e. negative and small, the neutron 
pair separation energies are expected to fall exponentially with increasing N, but the 
small value of T can bring this exponential fall very close to a linear one. 

The neutron pair separation energies of the even Sn isotopes from ^°^Sn to ^^°Sn 



(i.e. across the whole sdg neutron shell) have been fitted in E£3 using both theories. 
Furthermore, inila fit of the logarithms of the energies has been performed, since 
eq. (33.21) predicts a linear decrease of the logarithm of the energies with increasing 
N. Both fits give almost identical results. Eq. (33.21) (in which the free parameters 
are Gfl and T), gives a better result than eq. (33.20) (in which the free parameters are 
Vq and A/Vq) for every single isotope, without introducing any additional parameter, 
indicating that the higher order terms can be useful. 

One should, however, remark that ^^^Sn lies in the middle of the sdg neutron shell. 
Fitting the isotopes in the lower half of the shell (^°^Sn to ^^^Sn) and the isotopes in 
the upper half of the shell (^^^Sn to ^^°Sn) separately, one finds that both theories 
give indistinguishably good results in both regions. Therefore Q-deformation can be 
understood as expressing higher order correlations which manifest themselves in the 
form of particle-hole asymmetry. It is also known that a strong subshell closure exists 
at N=64 (which corresponds to ^^^Sn). The presence of this subshell closure can also 
affect the neutron pair separation energies in a way similar to the one shown by the 
data. 



In a a fit of the neutron pair separation energies of the Pb isotopes from ^'^^Pb to 
^°^Pb has also been attempted. In this case both theories give indistinguishably good 
fits. This result is in agreement with the Sn findings, since all of these Pb isotopes lie 
in the upper half of the pfh neutron shell. Unfortunately, no neutron pair separation 
energy data exist for Pb isotopes in the lower part of the pfh neutron shell. 



Concerning the values of T obtained in the case of the Sn isotopes (T = —0.0454, 
= —0.0447), one observes that they are shghtly smaller than the value (T = —0.0488) 
which would have been obtained by considering the neutrons up to the end of the sdg 
shell as lying in a single-j shell. This is, of course, a very gross approximation which 
should not be taken too seriously, since it ignores the fact that most properties of 
nuclei can be well accounted for by the valence nucleons alone, without being affected 
by the closed core. In the case of the Pb isotopes mentioned above, however, the best 
fit was obtained with T = —0.0276, which is again slightly smaller than the value of 
T = —0.0317 which corresponds to considering all the neutrons up to the end of the 
pfh shell as lying in a single-j shell. 

In summary, we have shown that pairing in a single-j shell can be described, up 
to first order corrections, by two Q-oscillators, one describing the J = pairs and the 
other corresponding to the J 7^ pairs, the deformation parameter T = InQ being 
related to the inverse of the size of the shell. These two oscillators can be used for 
forming an suq(2) algebra. A Hamiltonian giving the correct pairing energies up to 
first order corrections in the small parameter can be written in terms of the Casimir 
operators of the algebras appearing in the uq(2) D uq(1) chain, thus exhibiting a 
quantum algebraic dynamical symmetry. The additional terms introduced by the Q- 
oscillators serve in improving the description of the neutron pair separation energies 
of the Sn isotopes, with no extra parameter introduced. 

In the previous section a generalized deformed oscillator describing the correlated 
fermion pairs of J = exactly has been introduced. This generalized deformed 
oscillator is the same as the one giving the same spectrum as the Morse potential 
(sec. 35), up to a shift in the energy spectrum. The use of two generalized deformed 
oscillators for the description of J = pairs and J 7^ pairs in a way similar to the 
one of the present section is a straightforward task, while the construction out of them 
of a closed algebra analogous to the suq(2) obtained here is an open problem. The 
extension of the ideas presented here to the case of the BCS theory is an interesting 
open problem. 



33.3. Other approaches 



A g-deformed version of the pairing theory was assumed by Petrova £53 and Shelly 
Sharma with satisfactory results when compared to experimental data. The 



present construction offers some justification for this assumption, since in both cases 
the basic ingredient is the modification of eq. (33.4). It should be noticed, however, 
that the deformed version of eq. (33.4) considered in is different from the one 

obtained here (eq. (31.24)). A basic difference is that the deformed theory of Plr^^ 
reduces to the classical theory for g — > 1, so that g-deformation is introduced in order 
to describe additional correlations, while in the present formalism the Q-oscillators 
involved for Q ^ 1 reduce to usual harmonic oscillators, so that Q-deformation is 
introduced in order to attach to the oscillators the anharmonicity needed by the 
energy expression (eq. (33.6)). 



Continueing along the same line Shelly Sharma and Sharma c£3 derived Random 
Phase Approximation (RPA) equations for the pairing vibrations of nuclei differing by 
two nucleons in comparison to the initial one and applied their method to the study 
of the 0"^ states of the Pb isotopes, which offer a good example of pairing vibrations in 
nonsuper conducting nuclei. Furthermore, using deformed quasi-particle pairs coupled 
to zero angular momentum they developed a deformed version of the quasi-boson ap- 
proximation for 0"*" states in superconducting nuclei and tested it against a schematic 
two-level shell model B^. Another deformed two-level shell model has been developed 
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by Avancini and Menezes 



34. Anisotropic quantum harmonic oscillators with rational ratios of fre- 
quencies 

3-dim anisotropic harmonic oscillators with rational ratios of frequencies (RHOs) 
are of current interest because of their relevance as possible underlying syr nmetrj es of 
superdeformed and hyperdeformed nuclei E^'E^. In particular, it is thought I ^^^ J ^^^ that 
superdeformed nuclei correspond to a ratio of frequencies of 2:1, while hyperdef orrned 
nuclei correspond to a 3:1 ratio. In addition they have been recently connected P^lr^n 
to the underlying geometrical structure in the Bloch-Brink a-cluster model i^, and 
possibly to the interpretation of the observed shell structure in atomic clusters = 
especially after the realization that large deformations can occur in such systems 



The 2-dim RHO is also of interest, since its single particle level spectrum characterizes 
the underlying symmetry of "pancake" nuclei c^. 



RHOs are examples of maximally superintegrable systems P'^^r^l in N dimensions. 
Superintegrable systems in N dimensions have more than N independent integrals 
(constants of motion). Maximally superintegrable systems in N dimensions have 
2N— 1 independent integrals. 



anisotropic 



27^,279,28(1281 



The two-dim E^|269||27OH27i|,g73,g7|,0 three-dim 1300. 
harmonic oscillators have been the subject of several investigations, both at the clas- 
iicaLand the quantum mechanical level. The special cases with frequency ratios 1:2 
L23'E^ and 1:3 tH^I have also been considered. While at the classical level it is clear that 



the su(N) or sp(2N,R) algebras can be used for the description of the N-dimensional 
anisotropic oscillator, the situation at the quantum level, even in the two-dimensional 
case, is not as simple. 

In this section we are going to prove that a generalized deformed u(2) algebra is 
the symmetry algebra of the two-dimensional anisotropic quantum harmonic oscilla- 
tor, which is the oscillator describing the single-particle level spectrum of "pancake" 
nuclei, i.e. of triaxially deformed nuclei with to^ >> ^^y, i^z The method can be 
extended to the 3-dim RHO in a rather straightforward way. 



34-1- A deformed u(2) algebra 

Let us consider the system described by the Hamiltonian: 

H 



(34.1) 



where m and n are two natural numbers mutually prime ones, i.e. their great common 
divisor is gcd (m,n) = 1. 

We define the creation and annihilation operators 
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2' = ^ . 



x/m + ipx 
y/n + ipy 



V2 ' V2 
These operators satisfy the commutation relations: 



a, a 



Further defining 



1 

m' 
U -- 



1 

— , other commutators = 0. 

n 



W 



2 ^ ' J ' 2 
one can consider the enveloping algebra generated by the operators: 



6t 



So = -{U-W), H = U + W. 
These genarators satisfy the following relations: 

[5o,^±] = ±5±, [H,Si] = 0, for z = 0,±, 

and 



k=l \ 
m / 

s s^ = Y{ (u + 



2k - 1 



k=l 



2m 

2k -1' 
2m 



1=1 \ 

\{(w- 



21 



2n 

2e-V 

2n 



(34.2) 
(34.3) 

(34.4) 
(34.5) 

(34.6) 
(34.7) 

(34.8) 

(34.9) 
(34.10) 



The fact that the operators Si^ z = 0, ± are integrals of motion has been already 
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realized in 

The above relations mean that the harmonic oscillator of eq. (34.1) is described 
by the enveloping algebra of the generalization of the u(2) algebra formed by the 
generators S'o, 5"+, S_ and if, satisfying the commutation relations of eq. (34.8) and 



\S^, 5*+] — Fm^n{H, 5*0 + 1) — Fm,n{H, Sq 



(34.11) 



where 

So) = n [h/2 + So- n (^/2 ~So + . (34.12) 

In the case of m = 1, n = 1 this algebra is the usual u(2) algebra, and the operators 
5*0, S± satisfy the commutation relations of the ordinary u(2) algebra, since in this 
case one easily finds that 

[S_, S+] = -2So. (34.13) 

In the rest of the cases, the algebra is a deformed version of u(2), in which the 
commutator [S'-_,S'+] is a polynomial of 5*0 of order m + n — 1. In the case with 
m = 1, n = 2 one has 

5+] = 35o^ -HSo-^ + ^, (34.14) 

4 lo 

i.e. a polynomial quadratic in 5*0 occurs, while in the case of m = 1, n = 3 one finds 

5+] = -453 + - 1^0 - ^ + |, (34.15) 

i.e. a polynomial cubic in Sq is obtained. 

34-2. The representations 

The finite dimensional representation modules of this algebra can be found using 
the concept of the generalized deformed oscillator (sec. 12), in a method similar to the 



one used in liza'cra for the study of quantum superintegrable systems. The operators: 



= S+, A = S-, Af = So — u, u= constant, (34.16) 

where ti is a constant to be determined, are the generators of a deformed oscillator 
algebra: 

"Ar,^t] = ^t^ [^^^] = ^t^ = ^{HM), AA^ = <^{H,Af+l). (34.17) 

The structure function $ of this algebra is determined by the function Fm^n in eq. 
(34.12): 

^H,Ar)=Fm,n{H,Ar + u) = 

^ n {h/2 + H- " - — ) jn [m - - " + — ) . (34.18 

The deformed oscillator corresponding to the structure function of eq. (34.18) has an 
energy dependent Fock space of dimension + 1 if 

<^{E,0) = 0, ^{E,N +1) = 0, ^{E,k)>0, for A; = 1, 2, . . . , A^. (34.19) 



The Fock space is defined by: 



H\E,k>= E\E.k>, M\E,k>= k\E,k>, a\E,0>=0, (34.20) 

A''\E,k>^ ^J^E,k + l)\E,k + l>, A\E,k>^ ^J^E,k)\E,k-l > . (34.21) 
The basis of the Fock space is given by: 



,k>=^=(A^y\E,Q>, k^O,l,...N, (34.22) 



\E 



where the "factorial" [k] ! is defined by the recurrence relation: 

[0]! = 1, [A;]! = $(E, k)[k-l]\ . (34.23) 

Using the Fock basis we can find the matrix representation of the deformed oscillator 
and then the matrix representation of the algebra of eqs (34.8), (34.12). The solution 
of eqs (34.19) implies the following pairs of permitted values for the energy eigenvalue 
E and the constant u: 

where p = 1, 2, . . . , m, q = 1,2, . . . ,n, and 

2 V 2m 2n 7 ' ^ ^ 

the corresponding structure function being given by: 

^/r. N r ^ ^( 2p-l 2A;-1\"/.. 2g-l 2£-l\ 



j,^-^ y 2m 2m J iJ[\ 2n 2n J 

^ 1 r(ma:+p) T {{N - x)n + q + n) (34 26) 

m^n^ V {mx -\- p — m) T{{N — x)n + q) ' 

where T{x) denotes the usual Gamma-function. In all these equations one has N — 
0, 1, 2, . . ., while the dimensionality of the representation is given by A'"-|-l. Eq. (34.24) 
means that there are m • n energy eigenvalues corresponding to each N vahic. each 
eigenvalue having degeneracy + 1. (Later we shall see that the degenerate states 
corresponding to the same eigenvalue can be labelled by an "angular momentum".) 

It is useful to show at this point that a few special cases are in agreement with 
results already existing in the literature. 

i) In the case m = 1, n = 1 eq. (34.26) gives 

^E,x) = x{N + l-x), (34.27) 



while eq. (34.24) gives 



E = N + 1, 



(34.28) 



in agreement with Sec. IV. A of . 

ii) In the case m = 1, n = 2 one obtains for q = 2 



^E,x) =x{N + l-x) [N +l-x] , E = N + ^, (34.29) 



4' 



while for q = 1 one has 



<!>{E,x) =x{N + l-x){N+^-xY ^ = ^ + 1' (34.30) 



These are in agreement with the results obtained in Sec. IV.F of for the Holt 
potential (for (5 = 0). 

iii) In the case m = 1, n = 3 one has for g = 1 

<!>{E,x) =x{N +l~x) (^N + ^-x^ (^N +'^-xy E = N+'^, (34.31) 
while for g = 2 one obtains 

<^{E,x)=x{N+l-x)(^N+'^-x^(^N+^-xY E = N+ 1, (34.32) 
and for g = 3 one gets 

^{E,x) = x{N + 1- x) (^N x^ (^N x^ ^ = (34.33) 



These are in agreement with the results obtained in Sec. IV.D of for the Fokas- 
Lagerstrom potential. 

In all of the above cases we remark that the structure function has forms cor- 
responding to various versions of the generalized deformed parafermionic algebra of 
eq. (18.1), the relevant conditions of eq. (18.2) being satisfied in all cases. It is easy 
to see that the obtained algebra corresponds to this of the generalized parafermionic 
oscillator in all cases with frequency ratios 1 : n. 

The energy formula can be corroborated by using the corresponding Schrodinger 
equation. For the Hamiltonian of eq. (34.1) the eigenvalues of the Schrodinger equa- 
tion are given by: 

where n^. = 0, 1, . . . and n^^ = 0, 1, . . .. Comparing eqs (34.24) and (34.34) one con- 
cludes that: 

N = [n^/m] + [ny/n] , (34.35) 
where [x] is the integer part of the number x, and 

p = mod(n3;, m) + 1, q = mod{ny, n) + 1. (34.36) 



The eigenvectors of the Hamihonian can be parametrized by the dimensionahty 
of the representation A^, the numbers p, q, and the number k = 0,1, . . . , N. k can be 
identified as [n^/m]. One then has: 




34-3. The "angular momentum" quantum number 



It is worth noticing that the operators Sq, S± do not correspond to a generahza- 
tion of the angular momentum, 5*0 being the operator corresponding to the Fradkin 
operator Sxx — Syy F°°lr°^l. The corresponding "angular momentum" is defined by: 

Lo = {S+ - . (34.41) 

The "angular momentum" operator commutes with the Hamiltonian: 

[H,Lo]= 0. (34.42) 

Let \i > he the eigenvector of the operator Lq corresponding to the eigenvalue i. The 
general form of this eigenvector can be given by: 



\i>=J2 

k=0 



N 



(34.43) 



In order to find the eigenvalues of Lq and the coefficients we use the Lanczos 
algorithm i^, as formulated in From eqs (34.39) and (34.40) we find 



Lo\i>-- 



>-- 



' fc=0 



iv-1 t'^cj^f^jk + l) 



N 



k=l 



N 



(34.44) 



From this equation we find that: 



N 



Ck = {-lf2-^'^Hk{l/V2)/J\f, U^ = Y. 2""i^'(^/v^)/n! (34.45) 

n=0 

where the function Hk{x) is a generahzation of the "Hermite" polynomials (see also 
^■m), satisfying the recurrence relations: 

H^i{x) = 0, Ho{x) = 1, (34.46) 

Hk+iix) = 2xH,{x) - 2<!>l^^^{k)H,_,{x), (34.47) 
and the "angular momentum" eigenvalues i are the roots of the polynomial equation: 

HN+i{i/V2) = 0. (34.48) 

Therefore for a given value of N there are + 1 "angular momentum" eigenvalues 
i, symmetric around zero (i.e. if i is an "angular momentum" eigenvalue, then —i 
is also an "angular momentum" eigenvalue). In the case of the symmetric harmonic 
oscillator {m/n = 1/1) these eigenvalues are uniformly distributed and differ by 2. In 
the general case the "angular momentum" eigenvalues are non-uniformly distributed. 
For small values of analytical formulae for the "angular momentum" eigenvalues 
can be found E^. Remember that to each value of N correspond m ■ n energy levels, 
each with degeneracy + 1. 

In order to have a formalism corresponding to the one of the isotropic oscillator, 
let us introduce for every and (m, n, p, q) an ordering of the "angular momentum" 
eigenvalues 



nL,m,n,p,q 



where L = N and fx = -L, -L + 2, . . . , L - 2, L, (34.49) 



by assuming that: 

^L,m,n,p,g < ^L,m,n,p,g ^ ^ (34.50) 

the corresponding eigenstates being given by: 



^ Jl. (-2)'=fffc(£^'™™)/v^) 

\L,fi;m,n,p,q} = }_^ 



N 



L 
k=0 



^ k 



^0 M^J¥i^\ 

(34.51) 

The above vector elements constitute the analogue corresponding to the basis of 
"spherical harmonic" functions of the usual oscillator. The calculation of the "angular 
momentum" eigenvalues of eq. (34.49) and the coefficients di,d2, ■ ■ ■ yd^^i in the 
expansion of eq. (34.51) is a quite difficult task. The existence of general analytic 
expressions for these quantities is not obvious. The first few "angular momentum" 
eigenvalues are given by: 



nl,m,n,p,g _i_ 

'±1 — 



r(m + p) r(n + q) , , 

V ^i^j V ^y; (34.52) 



^ m'^n" r(p) r(g) 



and 

^0 



'±2 — ±1 



1 / r(m + p)r(2n + g) r(2m + p) r(n + g) \ 

^ ^m^n ^ r(p) r(n + g) r(m + p) r(g) / y ■ ) 

For L > 2 the analytic expressions of the angular momentum eigenvalues and the 
coefficients are longer, but their calculation is a straightforward task. Numerical 
results for these quantities in the cases of frequency ratios 1:2 and 1:3 are given in 
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After working out a few examples (see for details) one finds out the following 
points: 

i) In the basis described by eqs. (34.16)-(34.19) it is a trivial matter to distinguish 
the states belonging to the same irrep for any m : n ratio, while in the Cartesian basis 
this is true only in the 1:1 case. 

ii) In the 1:2 case the irreps have degeneracies 1, 1, 2, 2, 3, 3, 4, 4, ... , i.e. "two 
copies" of the u(2) degeneracies 1, 2, 3, 4, . . . are obtained. 

iii) In the 1:3 cases the degeneracies are 1, 1, 1, 2, 2, 2, 3, 3, 3, ... , i.e. "three 
copies" of the u(2) degeneracies are obtained. 

iv) It can be easily seen that the l:n case corresponds to "n copies" of the u(2) 
degeneracies. 

v) Cases with both m, n different from unity show more complicated degeneracy 
patterns, also correctly reproduced by the above formalism. In the 2:3 case, for 
example, the degeneracy pattern is 1, 1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 3, 3, ... . 

vi) The only requirement for each energy eigenvalue to correspond to one irrep of 
the algebra is that m and n have to be mutually prime numbers. If m and n possess a 
common divisor other than 1, then some energy eigenvalues will correspond to sums 
of irreps, i.e. to reducible representations. 

vii) The difference between the formalism used in i7^.i77[|78[|8i| ^^^^ ^j^g used 



here is that in the former case for given m and n appropriate operators have to be 
introduced separately for each set of {p, q) values, while in the present case only one 
set of operators is introduced. 



34.4^. Multisections of the isotropic oscillator 



In B2y the concept of bisection of an isotropic harmonic oscillator has been in- 
troduced. One can easily see that multisections (trisections, tetrasections, . . . ) can 
be introduced in a similar way. The degeneracies of the various anisotropic oscilla- 
tors can then be obtained from these of the isotropic oscillator by using appropriate 
multisections. 

Using the Cartesian notation {nx,ny) for the states of the isotropic harmonic 
oscillator we have the following list: 
N=0: (00) 
N=l: (10) (01) 



N=2: (20) (02) (11) 

N=3: (30) (03) (21) (12) 

N=4: (40) (04) (31) (13) (22) 

N=5: (50) (05) (41) (14) (32) (23), 
where = rix + Uy. The corresponding degeneracies are 1, 2, 3, 4, 5, 6, ... , i.e. these 
of u(2). 

A bisection can be made by choosing only the states with nj^=even. Then the 
following list is obtained: 





=0: 


(00) 


N= 


=1: 


(10) 


N= 


=2: 


(20) (02) 


N= 


=3: 


(30) (12) 


N= 


=4: 


(40) (04) (22) 


N= 


=5: 


(50) (14) (32) 



The degeneracies are 1, 1, 2, 2, 3, 3, . . . , i.e. these of the anisotropic oscillator 
with ratio of frequencies 1:2. The same degeneracies are obtained by choosing the 
states with ny=odd. Therefore a bisection of the isotropic oscillator, distinguishing 
states with mod(?T,y, 2) = and states with m^odijiy, 2) = 1, results in two interleaving 
sets of levels of the 1:2 oscillator. 

By analogy, a trisection can be made by distinguishing states with mod(nj^, 3) = 
0, or mod(ny, 3) = 1, or mod(ny, 3) = 2. One can easily see that in this case three 
interleaving sets of states of the 1:3 oscillator, having degeneracies 1, 1, 1, 2, 2, 2, 3, 
3, 3, ... , occur. 

Similarly a tetrasection can be made by distinguishing states with mod(nj^, 4) = 
0, or mod(nj^,4) = 1, or mod(ny,4) = 2, or mod(nj^,4) = 3. The result is four 
interleaving sets of states of the 1:4 oscillator, having degeneracies 1, 1, 1, 1, 2, 2, 2, 
2, 3, 3, 3, 3, .... 

By bisecting and trisecting Uy one is left with six interleaving sets of states 
with degeneracies 1, 1, 1, 1, 1, 2, 1, 2, 2, 2, 2, 3, 2, 3, 3, ... , i.e. degeneracies of the 
2:3 oscillator. 

By bisecting (or trisecting, tetresecting, etc) both and ny one is obtaining the 
original u(2) degeneracies of the isotropic oscillator. 

It is therefore clear that the degeneracies of all m : n oscillators can be obtained 
from these of the isotropic oscillator by appropriate multisections. In particular: 

i) The degeneracies of the 1 : n oscillator can be obtained from these of the 1:1 
(isotropic) oscillator by ra-secting Uy or n^. 

ii) The degeneracies of the m : n oscillator can be obtained from these of the 1:1 
oscillator by m-secting and ra-secting Uy. 



34-5. Connection to H^^^ 



1, n 



For the special case m 
received here coincides with the finite W algebra W 



2 it should be not iced t hat the deformed algebra 



(2) 
3 
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The commutation 



(2) 

relations of the W3 algebra are 

[Hw, Ew] = 2Ew, [Hw, Fw] = -2Fw, [Ew, F^] = + Cw, (34.55) 

[Cw, Ew] = [Cw. Fw] = [Cw. Hw] = 0, (34.56) 
while in the m = 1, n = 2 case one has the relations 

[AT, = A\ [AT, A] = -A, [A, A^ = SS'^ -^-HSo + ^, (34.57) 

4 16 

[H, A^ = [H, A] = [H, So] = 0, (34.58) 

with Sq = Af + u (where u a constant). It is easy to see that the two sets of commu- 
tation relations are equivalent by making the identifications 

Fw = aA\ Ew = pA, Hw = -2So + kH, Cw = f{H), (34.59) 

with 

P<y = \. k = \, f{H) = -^H' + \. (34.60) 



34-6. Discussion 

In conclusion, the two-dimensional anisotropic quantum harmonic oscillator with 
rational ratio of frequencies equal to m/n, is described dynamically by a deformed 
version of the u(2) Lie algebra, the order of this algebra being m + n — 1. The repre- 
sentation modules of this algebra can be generated by using the deformed oscillator 
algebra. The energy eigenvalues are calculated by the requirement of the existence of 
finite dimensional representation modules. An "angular momentum" operator useful 
for labelling degenerate states has also been constructed. The algebras obtained in the 
special cases with 1 : n ratios are shown to correspond to generalized parafermionic 
oscillators. In the special case of m : n = 1 : 2 the resulting algebra has been identi- 
fied as the finite W algebra W3 . Finally, it is demonstrated how the degeneracies of 
the various m : n oscillators can be obtained from these of the isotropic oscillator by 
appropriate multisections. 

The extension of the present method to the three-dimensional anisotropic quan- 
tum harmonic oscillator is already receiving attention, since it is of clear interest in 
the study of the_ayinmetries underlying the structure of superdeformed and hyperde- 



formed nuclei E^'E^ 



35. The use of quantum algebras in molecular structure 

The techniques developed in this article can be applied in very similar ways in 
describing properties of diatomic and polytomic molecules. A brief list will be given 
here. 



1) Rotatiqna 



sxiq{2) model 



108 



s nertra, 
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diatomic molecules have been described in terms of the 
I. As in the case of nuclei, g is a phase factor (g = e*"^). 



In molecules r is of the order of 0.01. The use of the sUq(2) symmetry leads to a 
partial summation of the Dunham expansion describing the rotational- vibrational 
spectra of diatomic molecules c£ j. M olecular backbending (bandcrossing) has also 
been described in this fra. niew or k Rotational spectra of symmetric top molecules 



have also been considered e22 
2 
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in the framework of the sUg(2) symmetry. 
Vibrational spectra of diatomic molecules have been described in terms 



g-deformed anharmonic oscillators having the suq(l,l) E23 or the Ug(2) D Og(2) 
symmetry, as well 



303 



3& 



„j J, .^^^^ terms of generalized deformed oscillators similar to the ones 

used in sec. 26 Sioliol T hese r esults, combined with 1), lead to the full summation 



of the Dunham expansion 



A two-parameter deformed anharmonic oscillator 
with Uqp(2) D Ogp(2) symmetry has also been considered 1^. 

3) The physical content of the anharmonic oscillators mentioned in 2) has been 
clarified by constructing WKB equivalent potentials (WKB-EPs) EIEhI and classical 
equivalent potentials @, similar to the ones of sec. 13. The results have been corrob- 
orated by the study of the relation between sug(l,l) and the anharmonic oscillator 



with x'^ anharminicities l23. The WKB-EP corresponding to the sug(l,l) anharmonic 
iscillator has been connected to a class of Quasi-Exactly Soluble Potentials (QESPs) 
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tial 



}_Generalized deformed oscillators giving the same spectrum as the Morse poten- 
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and the modified Poschl-' 



potential as well deformed oscillator 
have also been constructed. In addition, q- 



containing them as special cases 
deformed versions of the Morse potential have been given, either by using the sOg(2,l) 
symmetry tl3 or by solving a g-deformed Schrodinger equation for the usual Morse 
potential 



5) A g-deform ed ver sion of the vibron model for diatomic molecules has been 



constructed Eiil'liUEL^^^ in a way similar to that described in sec. 29. 

6) For vibrational spectra of polyatomic molecules a model of n coupled gener- 
alized deformed oscillators has been built i^^, containg the approach of lachello and 
Oss Eh^'EhI as a special case. 

7) Quasi-molecular resonances in the systems ^^C+^^C and ^^C-|-^^0 have been 
described in terms of a g-deformed oscillator plus a rigid rotator IhI. 

A review of several of the above topics, accompanied by a detailed and self- 
contained introduction to quantum algebras, has been given by Raychev 
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36. Outlook 



Nobody likes binding himself by statements concerning the future. However, we 
attempt to give here a partial list of open problems, roughly following the order of 
the material in this review: 

1 ) The list of physical systems which can be classified under a generalized deformed 
su(2) symmetry (sec. 17) or under a generalized deformed parafermionic oscillator 



scheme (sec. 18) can be enlarged. Self-similar potentials and isospectral oscillator 
Hamiltonian systems could probably be related to these symmetries. 

2) The description of B(E2) values in terms of the sUg(2) model attempted so far 
(sec. 21) takes into account only the kinematical deformation effects. In order to 
take into account dynamical deformation effects, one has to build a larger algebra, 
of which the quadrupole operators will be members. (This will then probably be a 
deformed version of an su(3) algebra.) These quadrupole operators should behave as 
irreducible tensors of rank 2 under sUg(2). 

3) The sug(2) prediction about B(E2) values increasing with increasing angular 
momentum J, supported by the predictions of other models as well (see sec. 21), 
requires further testing against detailed experimental data. 

4) The construction of Clebsch-Gordan coefficients for the subclass of generalized 
deformed su(2) algebras for which this could be possible (sees sees 17, 25) is an open 
problem. 

5) The sUg(3) D sOg(3) decomposition for sUg(3) irreps other than the completely 
symmetric ones (see subsec. 27.3 for the current state of the art) remains an open 
problem, the solution of which is necessary for developing a deformed version of the 
su(3) limit of the Interacting Boson Model. 

6) Realizations of multi-level shells models in terms of deformed bosons (see sees 
31-33 for some one- and two-level cases) should be further pursued. 

7) The symmetry algebras of the various 3-dim anisotropic harmonic oscillators 
with rational ratios of frequencies should be worked out, since they are of interest 
in relation to superdeformed and hyperdeformed nuclei, and possibly to deformed 
atomic clusters (see sec. 34 for references). A deformed u(3) algebra should occur in 
this case, which could serve as the basis for building a deformed analog of the Elliott 



model suitable for superdeformed nuclei. 
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In molecular physics (sec. 35) the study of vibrations of highly symmetric 
polyatomic molecules (including fullerenes) by these techniques is of interest. 
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